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FOREWORD 


This report is one of several to be published from research conducted 
under NASA Contract NAS8-26751 entitled "Reduction and Error Analysis of 
the AVE II Pilot Experiment Data. " This effort is sponsored by the NASA 
Office of Applications under the direction of Marshall Space Flight Center's 
Aerospace Environment Division. The results presented in this report 
represent only a portion of the total research effort. Other reports will 
be published as the research progresses. Previously published research re- 
sults appear in NASA TMX-64877 entitled "Data for NASA's AVE II Pilot 
Experiment, Part 1: 25-mb Sounding Data and Synoptic Charts." 
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REDUCTION AND ERROR ANALYSIS 
OF THE AVE II PILOT EXPERIMENT DATA 

by 

Henry E. Fuelberg* 

Center for Applied Geosciences 
Texas A&M University 

I . INTRODUCTION 

The pilot experiment for the second NASA Atmospheric Variability 
Experiment (AVE IIP) was conducted during a 24-hour period beginning 
at 1200 GMT on May 11, 1974, and ending at 1200 GMT on May 12, 1974. 
Fifty-four rawinsonde stations took part in the experiment supported 
by the National Aeronautics and Space Administration (NASA). Rawin- 
sonde soundings were taken at 3-hour intervals during the period 
which yielded a total of 470 soundings. The names and locations of 
the stations taking part in the AVE II pilot experiment are given in 
Fig. 1 and Table 1. 

This report describes the data reduction techniques used to process 
the original baseline, angle, and ordinate data, and presents an error 
analysis of the final output. Copies of all computer programs used in 
the reduction process are presented. Instructions for using the programs: 
and sample output are given as well. 


*Research Assistant, Center for Applied Geosciences, Texas A&M University 
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Rawinsonde 


Station Number 

11001 (MSF) 

22001 (OUN) 

22002 (FSI) 

22003 (LNS) 

22004 (FTC) 

22005 (CHIC) 

201 (EYW) 

202 (MIA) 
208 (CHS) 
211 (TPA) 
213 (AYS) 
221 (VPS) 
226 (MGM) 
232 (EVE) 
235 (JAN) 
240 (LCH) 
248 (SHV) 
250 (BRO) 
255 (VCT) 

260 (SEP) 

261 (DRT) 
265 (MAF) ' 
304 (HAT) 
311 (AHN) 
317 (GSO) 
327 (BNA) 
340 (LIT) 
349 (UMN) 
363 (AMA) 
402 (WAL) 
405 ( I AD ) 
425 (HTS) 
429 (DAY) 
433 (SLO) 
451 (DOC) 
456 (TOP) 
486 (JFK) 
494 (CHH) 

■ 518 (ALB) 

520 (PIT) 
528 (BUF) 
532 (PIA) 
553 (OMA) 
562 (LBF) 


Table 1 

Stations Participating in AVE II Pilot Experiment 


Location 

Marshall Space Flight Center, Alabama 

Norman, Oklahoma 

Ft. Sill, Oklahoma 

Lindsay, Oklahoma 

Ft. Cobb, Oklahoma 

Chickasha, Oklahoma 

Key West, Florida 

Miami, Florida 

Charleston, South Carolina 

Tampa, Florida 

Waycross, Georgia 

Eglin AFB, Florida 

Montgomery, Alabama 

Boothville, Louisiana 

Jackson, Mississippi 

Lake Charles, Louisiana 

Shreveport, Louisiana 

Brownsville, Texas 

Victoria, Texas 

Stephenville, Texas 

Del Rio, Texas 

Midland, Texas 

Hatteras, North Carolina 

Athe-ns, Georgia 

Greensboro, North Carolina 

Nashville, Tennessee 

Little Rock, Arkansas 

Monette, Missouri 

Amarillo, Texas 

Wallops Island, Virginia 

Dulles Airport, Virginia 

Huntington, West Virginia 

Dayton, Ohio 

Salem, Illinois 

Dodge City, Kansas 

Topeka, Kansas 

Kennedy Airport, New York 

Chatam, Massachusetts 

Albany, New York 

Pittsburg, Pennsylvania 

Buffalo, New York 

Peoria, Illinois 

Omaha, Nebraska 

North Platte, Nebraska 


Rawinsonde 


Station Number 

606 (PWM) 
637 (FNT) 
645 (GRB) 

654 (HUR) 

655 (STC) 
662 (RAP) 
712 (CAR) 
734 (SSM) 
747 (INL) 
764 (BIS) 


Table 1 (Continued) 

Stations Participating in AVE II Pilot Experiment 


Location 


Portland, Maine 

Flint, Michigan 

Green Bay, Wisconsin 

Huron, South Dakota 

St. Cloud, Minnesota 

Rapid City, South Dakota 

Caribou, Maine 

Sault Ste Marie, Michigan 

International Falls, Minnesota 

Bismarck, North Dakota 
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II. DATA REDUCTION PROCEDURES 

The steps used to reduce the angle and ordinate data to a finished, 
usable product are shown in Fig. 2. Each of these steps is described 
in the following sections of this report. General features of the 
programming language and the computer facilities used at Texas A&M 
University are described in Appendix A. 

A. Data Extraction 

The original records for all rawinsonde flights were forwarded to 
the Marshall Space Flight Center where angle, ordinate, and baseline 
data were extracted from the records and punched into cards by 
personnel from Texas ASM University under the supervision of NASA 
personnel. 

A baseline card was punched for each sounding which contained 
surface data and temperature and moisture calibration variables. 

Table 2 gives the name and column location of all data on the baseline 
cards. Figure 3 shows a sample baseline card. The station height was 
generally missing which was indicated by a field of nines in the 
appropriate columns. 

A series of data cards containing ordinate information was punched 
for each sounding. A typical sounding run ascended to about 20 mb and 
contained about 160 ordinate cards. A list of the data contained on 
an ordinate card is given in Table 3, and a sample card is shown in 
Fig. 3. Contact number and pressure were always indicated on the 
cards, but missing values of temperature ordinate,; humidity ord inace, 
or time vere indicated by filling the particular columns with nines', 
Information for the cards was obtained from the original strip chart 
of the sounding and the pres s\tre ; calibration chart for the individual 



Fig. 2. Steps of the data reduction process 
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Columns 

1 

3-6 

8-11 

13-18 

20-23 

25-28 

30-33 


35-38 

40-43 

45 - 4-7 

50-51 

53-54 

56-57 

59-62 

64-68 

72 


Table 2 


Baseline Data Cards 


Variable 

GMD type. 

Station elevation in meters. 

Pressure contact number at the, surface. 

Surface pressure in millibars. 

Surface temperature in degrees Celsius. 

Surface relative humidity in percent. 

Baseline temperature corresponding to an ordinate 
of 37.6 units. The temperature is always negative, 
but the minus sign is omitted. 

Baseline relative humidity at -40C and ordinate 46. 
Surface wind speed in meters per second. 

Surface wind direction in degrees. 

Month of rawinsondo release. 

Day of rawinsonde release. 

Year of rawinsonde release. 

Time of rawinsonde release in GMT. 

Station identification number. 

Angle identifier, A '1' indicates that angle data 
for the sounding run are given every minute while 
a blank indicates that angle data are available 
every 30 sec. This value was punched after the 
cards were received at Texas A&M University. 
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Co lumns 

1-3 

5-8 

10-13 

15-18 

20-24 

50-51 

53-54 

56-57 

59-62 

64-68 

72 

73 


Table 3 

Ordinate Data Cards 


Variable 

Pressure contact number corresponding to the given 
pressure . 

Pressure at the given contact. 

Tempera tu re o rd ina te . 

Humidity ordinate. 

Time from rawinsonde release. 

Month of rawinsonde release. 

Day of rawinsonde release. 

Year of rawinsonde release. 

Time of rawinsonde release in GMT. 

Station identification number. 

A ' 1 ' indicates the last card of the ordinate deck. 

A ' 1 ' indicates that either the temperature ordinate 
or the time from i-elease have been manually inter- 
polated. This value was punched after the cards 
were received at Texas A&M University. 
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sonde. If one wished to manually interpolate temperature ordinate or 
time, this could be done, and by punching a in column 73, the fact 
would be indicated by an asterisk in the final output. Manual inter- 
polation was rarely done, however. Interpolation schemes using the 
reduction program are described in detail in subsequent sections of the 
report. The date, time, and station number given in columns 50-68 are 
sometimes erroneous due to mistakes in keypunching and should not be 
relied upon by users, but this data is always correct on the baseline 
ca rd . 

Angle data for each sounding were punched on a separate deck of 
cards. Inmost cases, angles were available every 30 sec by reading 
them directly from the recorder charts, but in some instances, either 
at 0000 GMT or 1200 GMT, the rawinsonde operators recorded angles at 
only 1-min intervals. As shown in Table 2, this fact is indicated on 
the baseline card. Table 4 describes the contents of the angle cards 
while Fig. 3 shows a sample card. After investigation it was decided 
not to perforin manual smoothing or to manually check elevation angles, 
but instead to do this by use of the reduction program. Missing data 
were indicated by nines in the particular columns. As was the case for 
the ordinate data, the sounding descriptions given in columns 50-68 
are sometimes erroneous and should not be used. 

As angle and ordinate decks were completed the data were checked 
for errors by NASA personnel by computing first differences on the 
raw data. Errors were corrected as they were found; this was the 
first phase of error checking . 
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Columns 

1-5 

7-10 

12-16 

50-51 

53-54 

56-57 

59-62 

64-68 

72 

73 

74 


Table 4 

Angle Data Cards 


Variable 


Time from rawinsonde release. 

Elevation angle. 

Azimuth angle. 

Month of rawinsonde release. 

Day of rawinsonde release. 

Year of rawinsonde release. 

Time of rawinsonde release in GMT. 

Station identification number. 

A '1' indicates the last card of the angle dock. 

A '1' indicates that some type of manual smoothing 
was performed on the data at Texas A&M University. 
A '1' indicates that the elevation angle was less 
than some given quantity. 
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B . Transfer of die Data to Tape 

After the initial data check, the cards were forwarded to Texas 
A<SM University for further processing. A complete sounding deck 
consisting of a baseline card, ordinate cards, and angle cards was 
created for each rawinsonde run. Care was taken to insure that 
stations and times were correctly matched. The sounding decks were 
then arranged according to increasing station identification number, 
and within the same station by increasing time, in preparation for 
transfer to magnetic tape. 

A program was written to transfer the card decks to magnetic tape 
and to create a "leader card" which would precede the baseline informa- 
tion of each sounding. The "leader card" contains the date, time, 
and station number of the sounding that was found in columns 50-68 of 
the baseline card and also contains the total number of cards (records) 
for the sounding. This quantity enables one to quickly skip over a 
complete sounding when using the tape. The format of the leader card 
is given in Table 5; baseline, ordinate, and angle data on magnetic 
tape have the same format as when the data were on cards (Tables 2-4). 

A copy of the computer program and instructions for its use are found 
in Appendix B. 

A simple program was written to print the contents of the raw data 
tape. A copy of the program along with instructions for its use are 
found in Appendix C. Tlie program simply reads each record and prints 
the contents as they appeared on the original cards. The data list 
produced by the program proved valuable in later error checking 
procedures.;' 
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Co lumns 

1-2 

4-5 

7-8 

10-13 

15-19 

21-24 


Tab 1 5 

Format of the "^acler Cards" 


Variable 

Month of rawinsonde release. 

Day of rawinsonde release . 

Year of rawinsonde release. 

Time of rawinsonde release in GMT. 

Station identification number. 

Total number of cards (records) comprising the 
sounding. 
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C . Master Reduction Program 

This program was the most important in the data reduction process 
since it computed meteorological parameters for each pressure contact 
from baseline, ordinate, and angle data that had been previously stored 
on tape. A copy of the program is found in Appendix D along with 
instructions for its use. The accuracy of the output is described in 
Section III. The program can be conveniently divided into several 
functional components as follows to facilitate its description; a 
flowchart of the program is given in Fig. 4. 

1. Program Initialization . The first section of the program 
initializes the necessary arrays, defines functions to compute 
temperature (RESIST) and virtual temperature (XVIRT) , and defines 
constants to be used in humidity determination (DATAC) , All variables 
that will be printed and transferred to magnetic tape as final output 
are initially defined to be zero. The station roster containing the 
station identification number, elevation, latitude, longitude, and 
name is read from data cards followed by the number of soundings to 
process or skip. All of the remaining data is read from magnetic tape. 

2. Baseline Calibration , The first record to be read from tape 
is the "leader card" described in Table 5. If one desires to compute 
consecutive soundings beginning with the first sounding, the informa- 
tion on the "leader card" is not used; however, if one wishes to 
calculate certain soundings and skip others, the information is used 
to read raw sounding data, skip the computation phase, and then begin 
With the next "leader card" until the desired station and time are 
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located. The baseline data for the souiv. ag to be computed are read 
using the format described in Table 2. Baseline temperature constants 
are set up at this stage for later use in computing temperature at 
each pressure contact. The equations used in this section of the 
program and the temperature computation section are identical to those 
currently used by the National Weather Service (NWS) and are described 
by Billions (1965). The baseline equations are reproduced here: 


EXP jj. 6. 0082 991 - 0.9966256 • LQG e (2*37.6)] - 48000.0, (1) 


RTB 

RK1 = TO + 273.15, 

/ 

mi 


5.3018981 , 


RM2 

RM3 


~2.47991E-3+r(2.47991E-3) 2 -4.0-5.89986E-5-KMl]^ 
2.0-5. S9986E-5 


14000.0 

RTB 


EXP(RM2) 


( 2 ) 

(3) 

, and (4) 
(5) 


TO is the negative temperature corresponding to an ordinate of 37.6 
units; the value RM3 is used later. 

Constants to be used later in humidity computation are evaluated 
next. The scheme used in this program is based on a procedure 
described in IRIG Document 103-72 published by the White Sands Missile 
Range (1972). The 20 constants for humidity computation (HC1-HC20) 
are based on the following equations: 


HCl — Cl + C2 • RHO 

: : : : ( 6 ) 

HC20 = C39 + C40 • RHO , 
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Cl through C20 are constants determined by the type of humidity element 
used in the sonde, while RHO is the baseline humidity at a temperature 
of -40C and ordinate 46, The differences between the humidity 
computation procedures of the NWS and this program will be described 
in Section IXI-C. 

3. Reading the Ordinate Data , The ordinate data described in 
Table 3 are read from magnetic tape. Although pressure contact number 
and the corresponding pressure are always available, time from release, 
humidity ordinates, and temperature ordinates may be missing. Since 
time at each pressure contact is necessary in the reduction scheme, a 
linear interpolation procedure is vised to assign time where it is 
missing. If interpolation is performed, a locator array (IORIN) is 
used to mark this fact on the final output with an asterisk. Missing 
temperature values are interpolated in a later section, but humidity 

is not interpolated. 

4. Reading the Angle Data . Most raw angle data were obtained 
for 30-sec intervals, but at 0000 GMT and 1200 GMT, the data at some 
stations were available at only 1-min intervals (IlMIN = 1). In many 
cases the 30-sec angle data did not begin at the time of release but 
began at some time after release; therefore, the first angle record is 
checked to determine if it begins at time = 0. If it does not begin 
at time “0, a linear interpolation scheme is used to fill in the 
missing time* but the missing angles themselves are not interpolated. 
The fact that interpolation was performed would appear in the final 
output as missing wind data. Once the angles are determined to begin 
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Fig. 4. Flowchart of the Master Reduction Program. 
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at time of release, or after the missing times are interpolated back 
to the time of release, the angles are read without use of the 
procedure just described. This is possible because missing angle 
times are not encountered at any place in the data except at the 
beginning. The angles themselves may be missing, however. 

No angle data on the tapes that was at 1-min intervals failed to 
begin at the time of release; therefore, no interpolation was necessary 
for these data. The angle data are read at 1-min increments and then 
linearly interpolated to give values every 30 sec. This interpolation 
process is indicated by a statement at the top of the output data for 
that sounding. 

5. Assignment of Station Parameters . This section matches the 
station identification number of the sounding being computed with that 
of a station number in the roster so that station elevation and name 
are retrieved for further use. The subroutine MOPRT is called at 
this time to convert the numerical month of the sounding to the 
alphabetical name of the month. Also, azimuth angles for the four 
stations operated by the National Severe Storms Laboratory (22002-22005) 
are modified 180° to account for the different orientation of their 
rawinsonde tracking device. 

6. Temperature Computation . Temperature at each pressure contact 
is determined by use of the following equations (Billions, 1965): 


RE = EXP [l6. 0082991-0. 9966256‘LOG e (2*TORD)J - 48000.0, 


D/ T Art , RM3-RE , 
R4 14000.0 3 


1.0 


M “ 1.0/303.0 +4 . 6 7 7 4E - 4 • R4+1 . 1 1 2 7 8E - 5 • (R4)^ ’ and 


(7) 

( 8 ) 
(9) 
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TC = RK - 273.15 . 


( 10 ) 


TORD is the temperature ordinate at the pressure contact while RM3 is 
a result of the baseline computation (Eq. 5). TC is temperature in 
degrees Celsius. If the temperature ordinate is missing at a 
particular contact, the temperature is first defined to be missing and 
then linear interpolation is used in a subsequent section to fill the 
gap. The array IORIN is used to mark the location of this interpolation 
in the. final output. 

7. Computation of Other Thermodynamic Quantities . This section 
computes thermodynamic variables at each pressure contact; relative 
humidity is the first quantity to be computed. If the humidity ordinate 
is missing, all other moisture-related quantities at the contact are 
defined to be missing, and potential temperature (PTK) is computed 
according to the equation: 


PTK = TK 


1000.0 


0.286 


(ID 


where P is pressure and TK is temperature in degrees Kelvin at the 
contact. Virtual temperature is set equal to the ambient temperature 
when relative humidity cannot be computed. This equation and many of 
the other thermodynamic equations to follow may be found in texts such 
as Hess (1959) and TRIG Document 108-72 (1972). If the humidity 
ordinate is given, relative humidity (HP) is computed using the 
procedure described in IRIG Document 108-72 and given by: 


HP = (HC 1 -H 1 C 2 • TP) X 9 + (HC3+HC4 * TP) X 8 + ... + (HCl7-fHCl8*TP)x 

(HC19+HC20 • TP) . (12) 
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The values of TIC were obtained from the baseline calibration (Eq. 6); 
TP is temperature at the contact in degrees Celsius while. X is 
defined by the equation 


HORD-46 

41 


(13) 


HORD is the humidity ordinate at the pressure contact. If relative 
humidity is less than 57=, it and the other moisture variables are 
defined to be missing and Eq. (11) is used to compute potential 
temperature. In this situation, virtual temperature is again set equal 
to ambient air temperature. In some cases the procedure described 
above for computing relative humidity gives values greater than 1007.. 
This could be due to improper baseline temperature or humidity 
calibration or due to faulty sonde sensors. Minor errors in the 
baseline data can lead to the computed supersaturations . The reduction 
program truncates relative humidity to a maximum value of 100% for 
use in further computations, but prints the computed value even though 
it may be greater than. 1007,. Users may then decide what action to 
take regarding use of the data. 

The following equations are used to compute the remaining thermo- 
dynamic quantities based on relative humidity between 57, and 1007,. 

Vapor pressure in millibars (E) : 

E - HP ' • 0.0611 (10) ^ 7,5 ’ TP ^ 237,3+TP ^ (14) 


Mixing ratio in gm/Kg (W) : 
623 . 0 • E 

w — . . 

P-E 


( 15 ) 
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The specific heat of moist air in cal gm deg (CP): 


CP = 0.24(1.0 4- 0.84 


W 


1000.0 


) 


(16) 


The temperature of the dew point in degrees Celsius (TD) : 


TD = 


237.3-LO'G (E)-186.527 

e 

8.286-LOG (E) 
e 


(17) 


The potential temperature of moist air in degrees Kelvin (PTK) : 

RD/CP 


PTK = TIC 


1000. 0 

P-E 


(18) 


-1 -1 

RD is the gas constant for dry air, namely 6.87E-2 cal gm deg . 
Equivalent potential temperature in degrees Kelvin (EPOT) ; 


EPOT = PTK • EXP 


EL -W/ 1000.0 


(19) 


CP-TSA 

where EL is the Latent heat of vaporization at the lifting condensation 
level of air at the pressure contact, and TSA is the approximate 
temperature of air in degrees Kelvin at the lifting condensation level, 
of air at the pressure contact, and TSA is the approximate temperature 
of air in degrees JCe Ivin at the lifting condensation level. Virtual 
temperature in degrees Kelvin (TV); 


TV = 


TK 


1.0-(0.379*~) 


( 20 ) 


8. Computation of Height at 30-sec Intervals . Height of the sonde 
above ground level (YS) is computed using the hypsometric equation: 


YS = SYS + 


RD 


TV 4- TV , 
K K-l 


LOG 


'IC-1 


e\ P. 


K 


( 21 ) 
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where SYS is the height of the sonde at the preceding contact, g is 
acceleration of gravity, and K and K-l are the current and preceding 
pressure contacts, respectively. 

If there are less than ten angle observations for the sounding 
being computed, winds are not computed. Quantities dependent on the 
angles are then defined to be missing, and the program skips to the 
printing section of the program. 

If there are more than ten angle observations, height of the 
sonde is determined at 30-sec increments corresponding to the times 
of angle observations. This is achieved by a linear interpolation of 
height based on time from release which is common to both the height 
and angle data. 

9. Computation of Wind at 30-sec Intervals . Wind quantities are 

computed every 30 sec until the angle data are depleted. If angle 

data are not available at some time within the run, all quantities 

that are angle-dependent at that time are defined to be missing. The 

distance of the sonde from the station over a curved earth (SC) is 

computed using a procedure by Danielsen and Duquet (1966). The 

/ ' 

equation is : 

SC = RAD ^OS -1 

where RAD is the earth's radius, THETA is the elevation angle between 
the station and the sonde, and HINT is the height of the sonde above 
ground level. The X- and Y-location coordinates of the sonde, XS and 


COS (THETA) 

i hint 

1 RAD 


- THETA 


( 22 ) 
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ZS, respectively, are then easily determined from SC and the azimuth 
angle of the sonde (AZ) using the equations: 


XS = SC £siN(AZ)] , and 

(23) 

zs = sc £cos(az)] 

(24) 


Wind direction and speed at the surface are read from the baseline 
card; the east-west component (WWE ) and the north-south component (WSN) 
at the first level above the ground (30 sec after release) are obtained 
by computing a forward difference of the sonde locations over a 30-sec 
interval (DT = 30). After this time, centered differences over a 1-min 
interval (DT = 60) are tjged to obtain the components; in both cases 
the following equations are used: 


WWE(K-l) 


XS (K) -XS (K-2) 
DT 


, and 


(25) 


W sn(k-i) = gg.S L raster 


(26) 


The scalar wind speed is given by: 


WS (K-l) 



(27) 


Wind direction is obtained by first computing the value of A from the 
equation: 



WWE (K-l) 
WSN (K-l) 


57.29578 


(28) 


Since division by zero is undefined, steps must be taken to insure that 
this action is never attempted in Eq. (28). Therefore, the absolute 
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(29) 

(30) 

(31) 

(32) 

Finally, so wind direction (WD) is that direction from which the wind 
is blowing, the following corrections are applied: 

D = DIR +180°; (33) 

If D > 360°, WD = D - 360°; but ‘ (34) 

If D ^ 360°, WD = D. (35) 

10. Smoothing the Winds . A five point smoothing is performed on 
the wind components previously obtained at 30-sec intervals. The 
smoothing process does not begin until the time from release that 
corresponds to the sonde being 2.0 km above the surface. Once begun, 
the smoothing process continues through the remaining time periods 
except that if any one of the five values to be used in the smooth is 
missing, no smoothing is performed at that time, and the original 
components are retained. When smoothing is performed, the following 
equations ate used: 


value of WSN is never allowed to be smaller than 1.0E-5. After 
computing A, the following quadrant corrections are applied to 
determine wind direction: 


!f ~ £ 0 and WWE £ 0; 


DIR = 360° - A. 


WWK 

If £ 0 and WWE > 0; 
WSN 


If > 0 and WWE £ 0 ; 
WSN 


DIR = 180° - A. 
DIR = A + 180°. 


If —• > 0 and WWE > 0; 
WSN 


DIR « A. 
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Si = 

5 | 

E 

K-l 

WWE(K) • 

CD (kJ, and 

(36) 

S2 = 

5 

2 I 

K-l' 

[wSN(K) • 

CD(K)J. 

(37) 


Si is the smoothed east-west wind component, and S2 is the smoothed 
north-south wind component. The values of CD are the following 
binomial smoothing coefficients: 

CD(1) = 0.06, 

CD (2) - 0.25, 

CD (3) = 0.38, 

CD (4) ~ 0.25, and 
CD (5) - 0.06. 

The new components are used to determine a new wind direction and 
speed as was done in Eqs . (27-35). A locator array (I ASM) is used to 
indicate which values have been smoothed, and the array IEP is used 
to indicate which wind values have been based on elevation angles less 
than 9°. 

11. Interpolation of the Smoothed Winds to Pressure Contacts . 

Thus far the ordinate and wind data have been treated independently; 
the former were computed for pressure contacts which are approximately 
0.5 min apart near the surface and up to 1.5 min apart near the tops 
of typical soundings (20 mb). The latter were computed at 0.5-min 
intervals and then smoothed. Since time from release is common to both 
sets of datu, it is possible to interpolate the smoothed winds and 
assign values to times corresponding to the pressure contacts. The X 
and Y coordinates of the sonde are first linearly interpolated and used 
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to determine the range and azimuth angle of the sonde location. This 
is done in a manner similar to that of computing wind direction and 
speed except that coordinates and not wind components are used in this 
section. Wind components at the pressure contacts are determined by 
interpolation of the smoothed components which are then used to compute 
wind direction and speed at the contact according to Eqs. (27-35). If 
either of the two values used for linear interpolation is missing, the 
interpolated quantity is defined to be missing, and if either value 
has been smoothed or based on an elevation angle less than 9°, the 
interpolated quantity is defined to be smoothed and/or based on a low 
elevation angle, respectively. 

12. Printing the Results . This section of the program writes the 
results on paper and magnetic tape. Before doing so several last 
minute items are handled. Station elevation is added to the sonde 
heights above ground level which were used in wind computation. In 
the event that wind quantities became missing and could not be inter- 
polated before the last ordinate value, the winds at these levels are 
defined as missing. Arrays showing the location of temperature inter- 
polation, wind smoothing, and low elevation angles are defined to 
print an asterisk or blank instead of a 1 or a 0. 

Data for 45 pressure contacts are printed per page in addition to 
the headings. A description of the column headings is given in Table 6 
and a sample of the output is shown in Fig. 5. Station identification 
number, along with the date and time of release of the sonde are given 
at the top of the page. If angle data were available every minute 
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Table 6 

Column Headings for Output from the Master Reduction Program 


TIME (MIN) 
CNTCT 

HEIGHT (GPM) 

PRES (MB) 

TEMP (DG C) 

DEW PT (DG C) 
DIR (DG) 

SPEED (M./ SEC) 

U COMP (M/SEC) 

V COMP (M/SEC) 

POT T (DG K) 

E POT T (DG K) 
MX RTO (GM/KG) 
RH (PCT) 

RANGE (KM) 

AZ (DG) 


Time after balloon release. 

Contact number. 

Height of corresponding pressure surface in 
geopotential meters. 

Pressure in millibars. 

Ambient temperature in degrees Celsius. 

NOTE: An asterisk indicates that time from release 

and/or temperature were linearly interpolated. 

Dew point temperature in degrees Celsius. 

Wind direction measured clockwise from true north 
and is the direction from which the wind is b lowing. 

Scalar wind speed in meters per second. 

NOTE: An asterisk indicates that wind quantities 

are based on an elevation angle that is less than 9°, 

The E-W wind component, positive toward the east 
and negative toward the west. 

The N-S wind component, positive toward the north 
and negative toward the south. 

Potential temperature in degrees Kelvin. 

Equivalent potential temperature in degrees Kelvin, 

Mising ratio in grams per kilogram. 

Relative humidity in percent. 

Distance balloon is from release point along a 
radius vector. 

Direction toward balloon measured clockwise from 
true north. 


STATiCN NO. 522 
PEORIA, ILL 

It MAY 1974 

2315 GMT 156 26. 0 


TIME 

CNTCT 

HFIGHT 

PRES 

TEMP 

DEW PT 

DIR 

SPEED 

U COMP 

V COMP 

POT T 

E POT T 

MX RTO 

RH 

RANGE 

tl 

MIN 


GPM 

MB 

DC C 

OG C 

DG 

M/SEC 

M/SEC 

M/SEC 

DG K 

DG K 

GM/KG 

PCT 

KM 

OG 

0.0 ; 

6. 5 

200.0 

981.5 

19.4 

10.1 

290.0 

5.1 

4.8 

-1.7 

2 < ”5 . 2 

316.4 

8.0 

55.0 

0.0 

0. 

0.1 

7.0 

239.5 

977.0 

is. e 

9. 2 

308.0 

11.8 

9.3 

-7.3 

294.9 

314.9 

7.5 

53.7 

0.3 

114. 

0.4 

8. 0 

336. 7 

966.0 

18.6 

9. 3 

30 8.0 

11.8 

9.3 

-7.3 

295.6 

316.1 

7. 7 

54.8 

0.3 

114. 

0.8 

9.0 

434. 9 

95 5.0 

17. 5 

8.6 

306.6 

11.2 

9.0 

-6.7 

295.5 

315.2 

7.4 

55.7 

0.5 

120. 

1.0 

10.0 

524.7 

945 . 0 

16. S 

7.6 

305.5 

L 0. 7 

8.7 

-6.2 

295.3 

314.0 

7.0 

55.8 

0.6 

122. 

1.3 

11.0 

642.5 

932.0 

15.5 

6.7 

313.7 

10.9 

7.9 

-7.6 

295-3 

313.2 

6.6 

55.9 

0.8 

122. 

1.7 

12.0 

733.9 

922.0 

14. 5 

5.7 

316.8 

10.9 

7.5 

-7.9 

295.2 

312.1 

6.3 

55.8 

l.C 

126. 

2.1 

13.0 

826. 0 

912.0 

13.7 

6.9 

297.5 

10. 0 

8.8 

-4.6 

295.3 

313.8 

6.9 

63.6 

1.3 

130. 

2.3 

14.0 

918. 8 

902.0 

12.3 

5.2 

267.1 

11.3 

11.3 

0.6 

294.8 

311.4 

6. 2 

61.7 

1.4 

129. 

2.7 

15.0 

1012.2 

892.0 

11.7 

4.4 

246.9 

15.6 

14.3 

6.1 

295.1 

311.0 

5.9 

60.5 

1.6 

117. 

3.0 

16.0 

1125.4 

880.0 

10.8 

3.4 

246.3 

16.5 

15. 1 

6.6 

295.2 

310.4 

5.6 

60.2 

1.8 

105. 

3.4 

17.0 

1220.6 

870.0 

9.9 

2, 5 

265.7 

12.0 

12.0 

0.9 

295.2 

309.6 

5.3 

60.2 

2.1 

104. 

3.7 

18.0 

1316.6 

860.0 

8.9 

1. 6 

273.2 

11.3 

11.3 

-0.6 

295.1 

308.9 

5.0 

60.1 

2.3 

103. 

4.1 

19,0 

1413.4 

850.0 

7.8 

1.5 

274.8 

11.5 

11.4 

-1.0 

295.0 

308.7 

5.0 

64.1 

2.5 

102. 

4.4 

20.0 

1510.9 

840.0 

7.1 

0.9 

277.6 

11.5 

11.4 

-1.5 

295.2 

308.5 

4.9 

64.3 

2.8 

101. 

4.8 

21.0 

1629,1 

828.0 

5.7 

-0.4 

278.9 

11.7 

11.6 

-1.8 

294.9 

307.2 

4.5 

64.5 

3.1 

101. 

5.2 

22. 0 

1718. 5 

81 9 . 0 

5.0 

-0.3 

27 9,7 

11.7 

11.6 

-2.0 

295. 1 

30 1.7 

4.6 

68.2 

3.3 

101. 

5.4 

23.0 

1818.7 

809.0 

4.2 

-1.2 

280.3 

11.5 

11.4 

-2.1 

295.2 

307.2 

4.3 

68.0 

3.5 

101. 

5.B 

24.0 

1919.9 

799 1,0 

3.7 

o _ > 

231.8 

12.1 

11.9 

—2. 5 * 

295.5 

302.9 

2.6 

41.9 

3.8 

101. 

6.2 

25.0 

2012.0 

790,0 

4. 5 

-22.6 

280.9 

13.4 

13.2 

-2.5 

297.1 

299.5 

0.8 

11.8 

4. 1 

101. 

6.5 

26.0 

2136.5 

778,0 

4.5 

-22. 6 

278.8 

14.7 

14.6 

-2.3 

29 3.4 

300.9 

0.8 

11.8 

4.3 

101. 

6.8 

27.0 

2241.7 

768.0 

4.2 

-22.8 

276,1 

14.0 

13.9 

-1.5 

299.2 

301.6 

0.8 

11.8 

4.6 

101. 

7.2 

28.0 

2337.4 

739.0 

3.7 

-23.2 

272.1 

13.3 

13.3 

-0.5 

299.6 

302.0 

0.8 

11.9 

4.9 

100. 

7.5 

29.0 

2434. 1 

750.0 

3.3 

-23.4 

269.0 

13. 1 

13.1 

0.2 

300.3 

302.7 

0.8 

11.9 

5.1 

100. 

7.8 

30.0 

2 542 . 7 

740.0 

2.3 

-24. 1 

266.3 

13.2 

13.1 

0.8 

300.3 

302.6 

0. 7 

12.0 

5.4 

99. 

8.1 

31.0 

2630.4 

732.0 

2.0 

-24.3 

263.9 

13.4 

13.3 

1.4 

300.9 

303.2 

0.7 

12.0 

5.6 

99. 

8.5 

32.0 

2741.0 

722.0 

0.7 

-25.2 

261.7 

14.0 

13.8 

2.0 

300.6 

302.8 

0.7 

12.2 

5.9 

98. 

8.8 

33.0 

2841.6 

713.0 

0.2 

-25.6 

260.4 

14.7 

14.5 

2.4 

301.2 

303.3 

0.7 

12.2 

6.1 

97. 

9.1 

34.0 

2931,7 

705.0 

-1.1 

-26.5 

253.9 

15.4 

15. 1 

3.0 

300.7 

302.7 

0.6 

12.4 

6.4 

96. 

9.4 

35,0 

3034. 1 

696.0 

-1.6 

-26.8 

25 6.9 

16.0 

15.6 

3.6 

301.3 

303.2 

0.6 

12.4 

6.7 

96. 

9.7 

36. 0 

3160.6 

685.0 

-2.7 

-27.6 

254.5 

16.6 

16.0 

4.4 

301.4 

303.3 

0.6 

12.5 

7.0 

95. 

10.1 

37.0 

3265.3 

676.0 

-3.0 

-27.8 

251.4 

17.3 

16.4 

5.5 

302.2 

3 04.1 

0.6 

12.6 

7.3 

94. 

10.4 

38.0 

3371.2 

667.0 

-4.2 

-28.7 

249.4 

17.9 

16.8 

6.3 

302.0 

303.7 

0. 5 

12.7 

7.6 

93. 

10.7 

39.0 

3466.1 

£59.0 

-5.0 

-29.3 

247.3 

18.9 

17.4 

7.3 

302.2 

303.8 

0.5 

12.8 

7.9 

92. 

11.1 

40.0 

3574. 1 

650. 0 

-5.6 

-29.7 

244.3 

20. 5 

18.5 

8.9 

302.7 

304,3 

0. 5 

12.8 

8.3 

91. 

11.5 

41.0 

3695.5 

640.0 

-6. C 

-30.0 

240.8 

2 2.4 

19.6 

11.0 

303.5 

305.1 

0. 5 

12.9 

8.8 

89. 

11.8 

42.0 

3793.9 

632.0 

-6.5 

-27.8 

23 7,9 

23.3 

19.8 

12.4 

304.1 

306. 1 

0.6 

16.4 

9.2 

88. 

12.0 

43.0 

3905.8 

623.0 

-7.4 

-27.1 

236.2 

2 4.0 

19.9 

13.3 

304.3 

306.5 

0.7 

18.9 

9.5 

37. 

12.3 

44. 0 

4006. 5 

615.0 

-7. 2 

-24.1 

235.1 

24.2 

19.9 

13.9 

305.7 

308.5 

0.9 

24.5 

9.8 

85. 

12.9 

45.0 

4108.4 

607.0 

-7.7 

-22.1 

235.7 

25.3 

20.9 

14.3 

306.3 

309.7 

1.1 

30.2 

10.5 

83. 

13.2 

46.0 

4211,5 

599. 0 

-8.6 

-21.7 

236.9 

26.5 

22.2 

14.5 

306.4 

310.0 

1.1 

33.5 

10.9 

82. 

13.6 

47.0 

4315.6 

591.0 

-9.6 

-22.1 

238.1 

28.5 

24.2 

15. 1 

306.4 

309.9 

1.1 

34.9 

11.6 

81. 

12.9 

48.0 

4420.9 

533.0 

-9.4 

-23.4 

233.1 

30.1 

25.6 

15.9 

307. 8 

310.9 

1.0 

30.9 

12.1 

79. 

14.3 

49.0 

4527. 7 

575.0 

-9.4 

-25. 1 

236.8 

32.1 

26.8 

17.6 

309.0 

311.8 

0.9 

26.6 

12.8 

78. 

14.7 

50.0 

4635.9 

56 7.0 

-9.7 

-25.0 

234.3 

33.9 

27.5 

19.7 

309.9 

312.7 

0.9 

27.4 

13.5 

77. 
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Fig. 5. Sample output from the master reduction program 
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T|«E 

CNTCT 

HEIGHT 

PRES 

TEMP 

OEW PT 

or r 

SPEED 

U COMP 

V COMP 

POT T 

E POT T 

MX RTO 

RH 

RANGE 

AZ 

MIN 


GPM 

MB 

OG C 

DG C 

DG 

MY SEC 

M/SEC 

.M/SEC 

OG K 

OG K 

GM/KG 

PCT 

KM 

OG 

15.1 

5 L .0 

4773.1 

55 7.0 

- 10.3 

- 25.1 

231.4 

3 5-5 

27 . 8 

22.2 

310.8 

313.7 

0.9 

28.5 

14.3 

76 . 

15.4 

52.0 

4370.5 

550.0 

- 10.7 

- 22.9 

229.4 

36.5 

27.7 

23.3 

311.4 

314.9 

l.l 

35.8 

14.9 

75 . 

15 . 3 

53.0 

4982.9 

542.0 

- 11.7 

- 27.0 

227.5 

36.9 

27.2 

24.9 

311.5 

314.0 

0.8 

26.7 

15.7 

73 . 

1 6.2 

54.0 

5096. 9 

534.0 

- 11.4 

- 34.0 

226.4 

36.3 

26.3 

25.0 

313.1 

314.5 

0.4 

13.4 

16.5 

72 . 

16.4 

55.0 

5197.9 

527.0 

- 12.2 
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-4.0 

-2.8 

564.5 

999.9 

99.9 

999.9 

99.0 
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Fig . 5. ( Continued ) 
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(IlMIN = 1) instead of every one-half minute (IlMIN = 0), a statement 
to this effect is printed next. The three numbers in the upper 
right-hand corner are from left to right the number of pressure 
contacts printed, the pressure at the highest contact, and the value 
of IlMIN. Missing data are indicated by printing nines in the 
appropriate columns. The data transferred to magnetic tape are 
described in Table 7. 

D. First Differences of Contact Data 

A program was written to compute first differences of the output 
from the master reduction program in order to locate errors which had 
not been detected when first differences were computed on the raw data. 

A copy of the program is given in Appendix E. The program reads the 
sounding identification data and the contact data from magnetic tape 
and then computes first differences of each quantity except pressure 
contact number. If one of the values is missing, the first difference 
value is also defined to be missing; values marked with an asterisk in 
the contact data are also marked by an asterisk in the first difference 
program. Results from the first difference program are printed in the 
same format as the results from the master reduction program; a sample 
is shown in Fig. 6. Values printed at each contact are first 
differences except for pressure contact itself where the forward contact 
number is printed. The output from the program was scanned visually 
to locate possible errors. Actual errors were corrected by use of 
the following program. 


Table 7 
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Data Transferred to Magnetic Tape 


Quantity Dimension 

Station identification number None 

Station name 7 

Day of rawinsonde release None 

Month of rawinsonde release (IX, IY, IZ) None 

Year of rawinsonde release None 

Time of rawinsonde release None 

Number of pressure contacts in sounding None 

Minimum sounding pressure None 

Angle time identifier (IlMIN) None 

Time after balloon release 230 

Contact number 230 

Height of pressure surface 230 

Pressure 230 

Temperature 230 

Interpolation indicator 230 

Dew point temperature 230 

Wind direction 230 

Scalar wind speed 230 

Elevation angle indicator 230 

E-W wind component 230 

N-S wind component 230 

Potential temperature 230 

Equivalent potential temperature 230 

Mixing ratio 230 

Relative humidity 230 

Range of balloon 230 

Direction of balloon 230 
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TIME 

CNTCT 

HEIGHT 

PRES 

TEMP 

OEW PT 

MIN 


GPM 

MO 

DG C 

DG C 

0 . 1 

7.0 

39 . S 

- 4.5 

- 0.6 

- 0.9 

0.3 

8.0 

97.2 

- 11.0 

- 0 . 2 

0 . 1 

0.4 

9.0 

90.1 

- 11.0 

- 1.1 

- 0.7 

0.2 

io.o 

; 89.3 

- 10.0 

- 1.0 

- 1.0 

0.3 

11.0 

■ 117.8 

- 13.0 

- 1.0 

- 0 . 9 

0.4 

12.0 

91.4 

- 10. 0 

— 1,0 

- 1.0 

0.4 

13.0 

92 . 1 

- 10.0 

- 0 . 8 

1.2 

0.2 

14.0 

92.8 

- 10.0 

- 1.4 

- 1.7 

0.4 

15.0 

93.5 

- 10.0 

- 0 . 6 

- 0.8 

0.3 

16.0 

113.2 

- 12.0 

- 0.5 

- 1.0 

0.4 

17.0 

95.2 

- 10.0 

- 0.9 

- 0.9 

0,3 

13.0 

96 . 0 

- 10.0 

- 0.9 

"- 0.9 

0.4 

19.0 

96.7 

- 10.0 

- 1.1 

- 0 . 1 

0.3 

20.0 

97.6 

- 10.0 

- 0.7 

- 0.6 

0.4 

21.0 

118.2 

- 12.0 

- 1.4 

- 1.3 

0.4 

22.0 

89.4 

- 9.0 

- 0.7 

* 0 . 1 

0.2 

23 . 0 

100.2 

- 10.0 

- 0.9 

- 0.9 

0.4 

24.0 

101.2 

- 10.0 

- 0.5 

- 6.9 

0.4 

25.0 

92.1 

- 9.0 

0.5 

- 14.5 

0,3 

26 . 0 

124 . 5 

- 12,0 

0.0 

0.0 

0.3 

27.0 

105.2 

- 10.0 

- 0.3 

- 0.2 

0.4 

23.0 

95,7 

- 9.0 

- 0.5 

- 0.4 

0.3 

29.0 

96.7 

- 9.0 

- 0.3 

- 0.2 

0.3 

30.0 

108.6 

- 10.0 

- 1.0 

- 0 . 7 

0.3 

31 . 0 

87.7 

- 8.0 

- 0.3 

- 0.2 

0.4 

32.0 

110-6 

- 10.0 

- 1.3 

- 0.9 

0.3 

33.0 

100.6 

- 9.0 

- 0.5 

- 0.2 

0.3 

34,0 

90.2 

- 8.0 

- 1.3 

- 0.9 

0.3 

35.0 

102.3 

- 9.0 

- 0.5 

- 0.3 

0.3 

36.0 

126.5 

- 11.0 

- 1.1 

— 0 . R 

0.4 

37 . 0 

104 . 8 

- 9.0 

- 0.3 

- 0.2 

0.3 

38.0 

105.9 

- 9.0 

- 1.2 

- 0.9 

0.3 

39.0 

94.9 

- 8.0 

- 0.8 

- 0 . 5 

0.4 

40.0 

107,9 

- 9.0 

- 0.6 

- 0.4 

0.4 

41.0 

121. 5 

- 10.0 

- 0.5 

- 0.3 

0.3 

42.0 

98.4 

- 8.0 

- 0.4 

2.2 

0-2 

43,0 

111.9 

- 9.0 

- 0.9 

0.8 

0.3 

4 4.0 

100.7 

- 8,0 

0 . i 

3.0 

0.6 

45.0 

102.0 

- 3.0 

- 0.4 

1.9 

0.3 

46.0 

103.1 

- 8.0 

- 0.5 

0.4 

0.4 

47.0 

104.1 

- 8.0 

- 1.0 

- 0.4 

0.3 

48.0 

105 . 3 

- 8.0 

0 . 1 

- 1.3 

0.4 

49.0 

106 . 8 

- 8.0 

0.0 

- 1.7 

0.4 

50.0 

108.2 

- 8.0 

- 0.3 

0.1 

0.4 

51.0 

137.2 

- 10.0 

- 0.6 

- 0.1 


Fig. 6. Sample 
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DIP 

SPEED 

U COMP 

V COMP 

POT T 

E POT T 

MX PT " 

RH 

RANGE 

AZ 

DG 

M/SEC 

m/SEC 

M/SEC 

DG K 

DG K 

GM/KG 

per 

XM 

DG 

18.0 

6.7 

4.5 

- 5.5 

- 0.3 

- 1.5 

- 0.5 

- 1.3 

0.3 

114 , 

0.0 

0.0 

0.0 

0.0 

0.8 

1.2 

0 . 1 

1.1 

0.0 

0 . 

- 1.4 

- 0.6 

- 0.3 

0.6 

- 0.1 

- 0.9 

- 0.3 

0.9 

0.2 

6 

- 1.0 

- 0.4 

- 0.2 

0.4 

- 0.2 

- 1.2 

- 0.4 

0 , 1 

0.1 

2 . 

8.1 

0.2 

- 0.8 

- 1.3 

0.1 

- 0.8 

- 0.3 

0.1 

0.2 

0 . 

3.1 

- 0.1 

- 0.5 

- 0.4 

- 0.2 

- 1.1 

- 0.4 

- 0.1 

0.3 

4 

- 19.3 

- 0.9 

1.4 

3.3 

0.2 

1.7 

0.6 

7.8 

0.3 

4 . 

- 30.4 

1.4 

2.5 

5.2 

- 0.6 

- 2.4 

- 0.7 

- 2.0 

0.1 

- 2 . 

- 20.2 

4.2 

3.0 

5.5 

0.3 

- 0.4 

- 0.3 

- 1.1 

0.2 

-11 

- 0,6 

1.0 

0.8 

0.5 

0.1 

- 0.7 

- 0.3 

- 0.3 

0.2 

-12 

19.4 

- 4.5 

- 3.1 

- 5.7 

- 0.0 

- 0.8 

- 0.3 

0.0 

0.3 

-1 

7.5 

- 0.7 

- 0.7 

- 1.5 

- 0.0 

- 0.7 

- 0.3 

- 0.1 

0.2 

-t 

1.6 

0.2 

0.1 

- 0.3 

- 0.2 

- 0.1 

0.0 

4.0 

0.3 

-L 

2. .8 

- 0.0 

- 0 . 1 

- 0.6 

0.2 

- 0.2 

- 0.2 

0.2 

0.2 

-0 

1.3 

0.3 

0.2 

- 0.3 

- 0.3 

- 1.3 

- 0.4 

0.1 

0.3 

-0 

0.8 

- 0.0 

- 0.0 

- 0.2 

0.2 

0.4 

0.1 

3.8 

0.3 

-0 

0.5 

- 0.2 

- 0.2 

- 0.1 

0.1 

- 0.5 

- 0.2 

- 0.3 

0.2 

-0 

1.5 

0.6 

0.5 

- 0.4 

0.3 

- 4.3 

- 1.7 

- 26.1 

0.3 

0 

- 0.9 

1.3 

1.3 

- 0. 1 

1.7 

- 3.4 

- 1.8 

- 30.1 

0.3 

0 

- 2.1 

1.3 

1.4 

0.3 

1.3 

1.4 

o.o 

0.0 

0.3 

0 

- 2.8 

-o.e 

- 0.7 

O.B 

0.7 

0.7 

- 0.0 

0.0 

0.3 

-0 

- 4.0 

- 0.6 

- 0 . ( 

1.0 

0.5 

0.4 

- 0.0 

0-1 

0.3 

-0 

- 3.0 

- 0.2 

- 0.2 

0.7 

0.7 

0.6 

- 0 . 1 ) 

0.0 

0.2 

-0 

- 2.7 

0.0 

0.0 

0.6 

0.1 

- 0.1 

- 0,0 

0.1 

0.2 

-1 

- 2.4 

0.2 

0.1 

0.6 

0.6 

0.6 

- 0.0 

0.0 

0.2 

-1 

- 2.2 

0.6 

0.5 

0.6 

- 0.3 

- 0.4 

- 0.0 

0.1 

0.3 

-1 

- 1.3 

0.7 

0.7 

0.4 

0.5 

0,5 

- 0.0 

0.0 

0.3 

-1 

- 1.5 

0.7 

0.6 

0.5 

- 0.5 

- 0.6 

- 0.0 

0.1 

0.3 

1 

- 2.0 

0-7 

0.5 

0.7 

0.6 

0.5 

- 0.0 

0.0 

0.3 

~1 

- 2.3 

0.6 

0.4 

0.8 

0.1 

0.0 

- 0.0 

o.i 

0.3 

-1 

- 3.1 

Q .7 

0.4 

l.l 

0.8 

0.8 

- 0,0 

0.0 

0.4 

-1 

- 2.1 

0.6 

0.4 

0.8 

- 0.2 

- 0.3 

- 0.0 

0.1 

0.3 

-1 

- 2.1 

1.0 

0.6 

1.0 

0.2 

0. 1 

- 0.0 

0.1 

0.3 

-1 

- 2.9 

1,6 

1.1 

1.6 

0 . 5 

0.5 

- 0.0 

0 . 1 

0.4 

— 1 

- 3.5 

2.0 

l.l 

2,1 

0,3 

0.8 

- 0.0 

0.0 

0.4 

-1 

- 2.8 

0.9 

0.2 

1.4 

0.6 

1.0 

0. 1 

3.5 

0.4 

-1 

- 1.8 

0.6 

0,1 

1 .0 

0.2 

0.4 

0.1 

2.4 

0.3 

-1 

-l.l 

0.2 

- 0 . 1 

0.5 

1.3 

2.0 

0.2 

5.6 

0.4 

-2 

0.6 

1.1 

1.0 

0.4 

0.7 

1.2 

0.2 

5.7 

0.7 

-3 

1.2 

1.2 

1.3 

0.2 

0.1 

0,3 

0.1 

3.3 

0.4 

-1 

1.2 

2.0 

2.0 

0.6 

- 0.0 

- 0.1 

' - 0.0 

1.5 

0.6 

-1 

0.1 

1.6 

1.4 

0.8 

L . 4 

1 . 1 

- 0. 1 

- 4.0 

0.5 

-l 

- 1.4 

1.9 

1.2 

1.7 

1.2 

0.8 

- 0.1 

- 4.3 

0.7 

-1 

- 2 . 4 ' 

1,8 

0.7 

2.2 

0.9 

1.0 

0.0 

0.8 

0.7 

-1 

- 2.9 

1.7 

0.3 

2.4 

0.9 

0.9 

0.0 

1.1 

0.3 

-1 


tput from. the first difference program. 
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time 

CNTCT 

HEIGHT 

PRES 

TEMP 

OEw P’ 

Dt*- 

SPEEO 

a ccmp 

V C IMP 

POT T 

E PCT T 

MX GTG 

RH 

RANGE 

it 

MIN 


GPM 

MB: ' 

DG C 

OG C 

G5 

M/SEC 

M/SEC 

M/SEC 

OG * 

OG K 

GM/KG 

PCT 

KM 

OG 

0.3 

52.0 

97.3 

-7.0 

-0.4 

2.1 

-2.0 

1.0 

-0.1 

1.6 

0.6 

1.3 

0.2 

7.3 

0.6 

-l. 

0.4 

53.0 

112.5 

-3.0 

-1.0 

-4.1 

- 1.8 

0.4 

- 0.5 

t . 1 

O.l 

-0.9 

-0.3 

-9.1 

0.8 

-2. 

0.5 

54.0 

114.0 

-8.0 

0.3 

-6.9 

-I.l 

- 0.6 

- 1.0 

O.l 

1.6 

0.5 

-0.4 

- 13.3 

0.3 

-1. 

0.2 

55.0 

101.0 

-7.0 

-0.7 

-0.5 

-0.4 

-0.7 

-0.6 

- 0.3 

0.3 

0.3 

- 0.0 

O.l 

0.5 

-I. 

0.5 

5o.O 

1-31.5 

-9.0 

-1.0 

-0.7 

0.1 

-1.9 

-1.4 

-1.4 

0.4 

0.3 

- 0.0 

0. 1 

0.9 

-1. 

0.5 

57.0 

lia.4 

- 6.0 

- 0.8 

- 0.6 

l.l 

-Q.5 

0.1 

-0.3 

0.4 

0.3 

- 0.0 

0.1 

0.7 

-1. 

0.4 

53.0 

104.9 

-7.0 

0.0 

0.0 

1.5 

0.3 

o. a 

-0.5 

L.2 

1.3 

0.0 

0.0 

0.7 

-1. 

0.3 

59.0 

106.3 

-7.0 

-0.6 

-0.4 

1.2 

0.5 

0.8 

-0.2 

0.6 

0.6 

-0.0 

0.1 

0.6 

-1. 

0.4 

60.0 

123. 1 

- 8.0 

-0.6 

-0.4 

1.3 

-0.4 

0-1 

-0.9 

0.8 

0.8 

-0.0 

0.1 

0.9 

-1. 

0.4 

61.0 

109.0 

-7.0 

-l.t 

- 0.8 

1.2 

-0.7 

-o.l 

-1.0 

-0.1 

-0.1 

-0.0 

O.l 

0.7 

-1. 

0.3 

62.0 

110 . 1 

-7.0 

-I. 1 

- 0.8 

0.9 

-0.3 

0.1 

-0.6 

-0.0 

- 0 . 1 

-0.0 

0.1 

0.5 

-0. 

0.4 

63. C 

111.4 

-7.0 

-0.7 

-0.5 

1.1 

0.6 

0.9 

-O.l 

0.5 

0.5 

-0.0 

O.l 

0.7 

-0. 

0.4 

64.0 

112.7 

-7.0 

-0.9 

-0.7 

0.9 

2.1 

2.0 

0.3 

0.2 

0.1 

-0.0 

0.1 

0.3 

-0. 

0.4 

65.0 

97.6 

-6.0 

-1.1 

- 0.8 - 

0.7 

1.9 

1.8 

0.7 

-0.2 

-0.2 

-0.0 

O.l 

0.9 

-a. 

0.4 

6o . 0 

146.2 

-9.0 

-1.1 

- 0.8 

0.7 

1.2 

1.2 

0.3 

0.5 

0.4 

- 0.0 

O.l 

0.9 

-0. 

0.4 

67.0 

ii6.8 

-7.0 

-1.2 

-0*9 

0.7 

0.3 

0.5 

-0.2 

-O.l 

-O.l 

-0.0 

0.1 

1.0 

- 0 . 

0.3 

63.0 

118.2 

-7.0 

-0.9 

-0.7 

0.7 

-0.3 

-0.0 

-0.6 

0.3 

0.2 

-0.0 

0.1 

0. 7 

- 0 . 

0.3 

69.0 

102.5 

- 6.0 

-0.9 

-0.7 

0.7 

-0.3 

0.0 

-0.5 

0.1 

0.0 

-0.0 

0.1 

0.7 

-0. 

0.4 

70. 0 

120.9 

-7.0 

-1.2 

-0.9 

0.6 

-0.2 

0.0 

-0.5 

-0.0 

-O.l 

-0.0 

0.1 

0.9 

-0. 

0.4 

71.0 

140.1 

- 8.0 

-0.7 

-0.5 

-0.1 

0.1 

0.0 

0.1 

0.9 

0.9 

-0.0 

O.l 

0.9 

-0. 

0.4 

72.0 

124.5 

-7.0 

-0.8 

-0.6 

-l.l 

0.3 

-o.l 

0.8 

0.5 

0.5 

-0.0 

0.1 

0.9 

-0. 

0.4 

73.0 

106.1 

-6.0 

-0.8 

-0.6 

-1.2 

0.3 

0.2 

1.1 

0.3 

0.3 

-0.0 

0.1 

0.9 

-0. 

0.4 

74.0 

127.8 

-7.0 

-0. 5 

-0.4 

-1.0 

1.5 

0.9 

1.4 

0.9 

0.9 

-0.0 

O.l 

0.9 

-0. 

0.4 

75.0 

111.1 

- 6.0 
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E. Correct Errors on Tape 

Many errors in the baseline, angle, and ordinate data were found 

after output from the first difference program was examined. Since 

most soundings contained at least one error, a program was written to 

substitute correct data records for incorrect ones. A copy of the 

program is presented in Appendix F. The original data were read from 

tape; the corrections were read from cards; and the corrected data 

were then transferred to a second tape. Three types of cards were 

needed to make the corrections. The first card specified the number 

of records to be corrected within the given sounding. If no corrections 

were needed, the data were simply transferred from the first tape to 

* ■ 

the second tape without modification. If corrections were needed, a 
second card was read for each correction that specified the location 
of the record to be corrected. This location was determined by manually 
counting down to the error on the raw data listing. For example, an 
error may be located on the fifth record of the sounding. Finally, a 
corrected data card with the original format followed each error 
location card. If corrections were needed in a sounding, the program 
read each record from the original tape and transferred it to the new 
tape until the location corresponding to the correction was reached; 
then the corrected card was substituted for the old record, and this 
was transferred to tape. The process continued until all corrections 
to the sounding had been made and all records of the original tape had 
been transferred to the new tape. 


F. Master Reduction Program Using Corrected Data 

Since many corrections were made to the original data, it was 
decided to recompute each sounding using the master reduction program 
described in Section II-C. No changes were made in the program or the 
type of output produced. 

G . Interpolation of Pressure Contact Data to 25-mb Increments 

Meteorological charts are frequently plotted for constant pressure 
surfaces; therefore, a program was written to interpolate the output 
from the master reduction program onto surfaces that are 25 mb apart. 

A copy of the program is found in Appendix G, and a flowchart is 
given in Fig. 7. The program assigns all of the surface values to 
level one (I = 1) of the new arrays (B and IB) and then linearly 
interpolates quantities to 25-mb pressure surfaces from 1000 mb (I = 2) 
up to 25 mb (I = 41). Since pressure contacts used in. interpolation 
are only about 12 mb apart near the surface and 3. mb apart near the 
top of the sounding, it is possible to use linear interpolation 
instead of logarithmic ' interpolation. If the particular 25 -mb level 
does not lie between the surface pressure value and the minimum pressure 
at the top of the sounding, the data for that 25-mb pressure surface 
are defined to be missing. Since wind direction may oscillate around 
360°, care must be taken to guard against fictitious directions. If 
data needed in the interpolation of a quantity are missing, the inter- 
polated quantity is also defined to be missing, and if quantities marked 
by an asterisk are used in the interpolation process, they are also 
marked by an asterisk in the output of this program, A sample of the 




























Fig. 7. (Continued) 
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output is presented in Fig. 8 where the column headings are the same 
as those described in Table 6. The results of the interpolation 
program are transferred to magnetic tape in the same order as described 
in Table 7, except items which had a dimension of 230 now have a 
dimension of 50. 

H. Errors Remaining in the Output 

In spite of efforts to correct all errors in the final output, 
some, remaining errors have been detected. Since the errors are few 
in number, and since another pass through the data would be a lengthy 
and costly process, these errors were not changed. They are listed 
in Appendix H along with the necessary steps to correct each of them 
where possible. Missing soundings are listed in Appendix I. 
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Fig. 8. Sample output from the program to convert pressure contact data to 


25 -mb increments . 
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III. ERROR ANALYSIS OF THE DATA 

Accuracy of meteorological data is an important consideration in 
its use. This is especially true when small-scale temporal and spatial 
features are being studied since these features may have amplitudes 
which approach the limit of accuracy of the data. The accuracy of 
data is dependent on several factors -- the type of equipment being 
used at a station, the reduction procedures used to process the data, 
human factors in calibrating and tracking the sonde, and handling the 
data. Current estimates of errors in the input data will be given in 
this section along with a discussion of how these errors affect the 
final product of the master reduction program and the program to 
interpolate data to constant-pressure surfaces. 

A. Stated Accuracy of the Thermodynamic Data 

1. Temperature . Case (1962) adopts a root-mean-square (rms) error 
of 0.7C based on laboratory determination, while Weidner and Chambers 
(1967) believe that a value of 1.4C is more realistic operationally. 
Hodge and Harmantas (1965) found rms differences between military and 
NWS sondes to be about 0.5C under field conditions. It is believed that 
rms errors in the AVE IIP temperature data are approximately 1C. 

2. Pressure . Radiosondes are calibrated after manufacture so that 
each pressure contact on the baroswitch corresponds to a certain 
atmospheric pressure. By comparing baroswitch pressures with those 
obtained hydrostatically using radiosonde temperatures and balloon 
heights measured by a high precision radar, Weidner and Chambers (1967) 
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found that baro switch errors are generally small. The rms errors 
were determined to be approximately 1.3 mb from the surface to 400 mb, 

1.1 mb from 400 to 100 mb, and 0.7 mb between 100 and 10 mb. These 
errors are considered representative of the AVE IIP data. 

3. Humidity . A carbon strip element is used now as the humidity 
sensor in radiosondes. Case (1962) reports an rms error of 10% in the 
sensor. The current specifications for the carbon element are +37° at 
25 C, +5% at temperature less than 25C plus an allowance for hysteresis 

of j4% above a value 33% relative humidity, and +5% at humidities between 
10% and 337o (Brousaides, 1973). Brousaides concludes that use of 
carbon humidity data at values below 20% is unwarranted except for the 
trend information that is provided. AVE IIP humidity data are thought 
to be generally within 10% of the true values. 

4. Pressure Altitude . Altitude of the sonde is determined by 
use of the hypsometric equation in which pressure, temperature, and 
humidity are variables. Errors in temperature, however, are the 
greatest source of error in height determination. The worst possible 

situation arises when temperature errors are systematically of the same 

* 

sign through great depths of the atmosphere. Table 8 shows the height 
error at various pressures due to a systematic temperature error of 1C. 
These errors are larger than would be expected in the AVE IIP data 
inasmuch as the errors in temperature generally are of both signs. 
Scoggins and Smith (1973), based on a study of AVE I data, concluded 
that realistic rms errors in pressure altitude are on the order of 10 
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Table 8 


Errors in Altitude Due to a Systematic Temperature 

Error of 1C 


Pressure 

Altitude 

mb 

SPH 

50 

87.7 

100 

67.5 

200 

47.1 

300 

35.3 

400 

26.8 

500 

20.7 

600 

15.3 

700. 

10.7 

800 

6.7 

900 

3.1 

1000 

0 
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gpm at 500 mb, 20 gpm at 300 mb, and 50 gpm at 50 mb. An analysis of 
the AVE IIP data suggests comparable errors. 


B. Accuracy of the Wind Data 

The AVE IIP wind data were computed in several stages beginning 
with finite differencing over 1-min intervals each 30 sec. These 
"raw" winds were smoothed, and the smoothed values were then inter- 
polated to pressure contacts. The accuracy of the "raw" wind data will 
be described in detail, but the effects of smoothing and interpolation 
on accuracy must be described in more general terms. 

1. Accuracy of the "Raw" Winds . Scoggins (1963) developed a 
statistical technique for evaluating errors in wind data that will be 
used in this section. If F is a function of X, Y, and Z, and if these 
are in error by AX, AY, and AZ, then F is in error by an amount AF 
which is given by the Taylor's series as: 


AF 


dF 

ax 


AX + 


oF 

dY 


AY + 


dF 

dZ 


AZ 


(38) 


The equation is valid if AX, AY, and Az are small compared with Sf/Sx, 
dF/oY. and dF/Sz; it was assumed that all terms involving squares, 
higher powers, and cross products of AX, AY, and AZ are negligible, 
and that higher derivatives are small in value. After squaring both 
sides of Eq. (38) and replacing each term by its average value, the 
following equation is obtained when one assumes that errors in X, Y, 
and Z are independent: 


2 .dF .2 , .dF ,2 , .dF N 2 

^dx °x^ W ct y^ + ^dz °zV 

The symbol cr denotes rms error. 


(39) 


Radar or radio direction- finding equipment is used to measure 
azimuth and elevation angles of the balloon while the hypsometric 
equation is used to determine its height. For the purpose of this 
discussion, a flat earth will be assumed although this was not done in 
the actual reduction process described earlier. The X, Y, and Z 
coordinates of the balloon are given by: 

X - Z COT(E) SIN(A) (40) 

Y = Z COT(E) COS (A) (41) 

Z = - T* LOG ) (42) 

g e\p 2 / 

where E is elevation angle, A is azimuth angle, R is the gas constant, 
g is acceleration due to gravity, T* is mean virtual temperature, and 
p^ and p 2 are pressures at two levels. Errors in X, Y, and Z can then 
be evaluated by use of Eq. (39). The assumptions made are that errors 
in E, A, T *, and p are independent of each other, and that errors in 
each parameter are normally distributed. Results from the AVE I angle 
data, recorded at 6-sec intervals, indicate that there is some systematic 
component to the elevation angles at values below 10° , but generally 

the above assumptions are thought to be valid for the AVE IIP data 

which was recorded at 30-sec intervals. Errors in Z are mostly due to 
errors in temperature and not errors in pressure or moisture. Applica- 
tion of Eq. (39) to Eqs. (40-42) yields: 

’•Personal communication with Mr. W. W. Vaughan, Chief, Aerospace 
Environment Division, NASA Marshall Space Flight Center, Huntsville, 
Alabama, 35812. 
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( 43 ) 


(44) 


(cr S L o G ~i ) 2 
T g e p ' 


(45) 


It is generally assumed that cr, = cr - 0.05° for the GMD-1 used in the 
' A E 

AVE II pilot experiment (Danielsen and Duquet, 1966); however, Weiss 

and Georgian (1969) indicate that these values are not equal, are 

larger than those generally assumed, and greatly dependent on elevation 

angle. A value of 1C is assumed for If one assumes that A = 90°, 

the value of cr will be a maximum, while if A = 0°, the value of cr is 
A Y 

a maximum. The rms error of wind speed components cr and a is given 

wx wy 


by 


a 


wx 


a 


fa 2 -ter 2 ]* 
L x2 x 1 _ I 

At 


" , and 


wy 


(46) 

(47) 


where subscripts 2 and 1 refer to consecutive measurements. These 

2 2 2 2 

equations can be simplified by assuming a „ = a and cr = cr which 

xx Y y x 

is reasonable over short time periods. The "raw" winds computed by 
the AVE IIP reduction program are based on angle measurements that are 
60 sec apart (At = 60). Finally, the rms error of the scalar wind (cr ) 


is given by: 
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<? = T^ 2 + a 2 1 ^ . (48) 

Ivl |_ wx wyj 

a , is a maximum at either 0° or 90°. 

Ivl 

Figure 9 shows the results of the computations at various eleva- 
tion angles and pressures; heights used for corresponding pressure 
Values are those of a standard atmosphere. RMS errors increase with 
increasing height and with decreasing elevation angle. Maximum rms 
errors for AVE IIP data at 700 mb range from about 2.5 m/sec at an 
elevation angle of 10° to about 0.5 m/sec at an elevation angle of 40°. 
At 500 mb the errors are 4.5 m/sec and 0.8 m/sec at the same elevation 

angles; while at 300 mb the errors are 7.8 m/sec and 1.0 m/sec, 

respectively. Users of the data should remember that these values are 
maxima since the azimuth angle was assumed to be either 0" or 90°. 

These errors of the "raw" wind are in agreement with those cited 
for the AVE I data by Scoggins and Smith (1973 ). Data published by 
the Air Force Missile Test Center (1963) are also in good agreement 
with those of the AVE II pilot experiment. RMS vector errors were given 

as a function of altitude and constant mean wind speed through a layer. 

By assuming an ascent rate of 1,000 ft/min, the elevation angle of the 
sonde at each altitude could be determined. RMS vector errors for an 
elevation angle of 12.5° are plotted on Fig. 9. The results agree well 
with the. AVE II results. Results by Reiter (1963), given in Fig. 9 
for an elevation angle of 10°, are also consistent with those of the 
AVK [I pilot experiment. The errors cited by Reiter are based on a time 
interval of 2 min and an error in angle measurement of 0.075°. 



When 8 is small, sin 8 # tan 8; therefore, rms errors in wind 


direction (cr dir) may be approximated by: 



where Ivl is the scalar wind speed. Since the error is .dependent on 
the azimuth angle of the sonde, maximum errors can be given by making 
A = 0°. Results for the AVE IIP data are plotted as a function of 
pressure and elevation angle in Fig. 10. At 700 mb the scalar wind 
speed is assumed to be 15 m/sec, while at 500 mb, 300 mb, and 100 mb, 
it is assumed to be 20 m/sec, 25 m/sec, and 20 m/sec, respectively. 

Errors for other geometric conditions and wind speeds would be different. 
The errors increase with decreasing pressure and decreasing elevation 
angle. Errors in AVE IIP wind direction at 700 mb range from about 
9.5° at an elevation angle of 10° to about 1.3° at an elevation angle 
of 40°. At 500 mb the errors are 13,4° and 1.8 at the same elevation 
angle, while at 300 mb the errors are 18.0 and 2.5, respectively. If 
the azimuth angle at 300 mb were 90° instead of 0°, the rms error would 
be 2.4° at an elevation angle of 10°. The results are in agreement 
With those cited for the AVE I data by Scoggins and Smith (1973), and 
those obtained ay personal communication with Mr. W. W. Vaughan of 
NASA - Huntsville, Alabama. If one assumes that the rms vector errors 
given by the Air Force Missile Test Center (1963) are maximum errors, 
the maximum rms direction error will be given by: 

= TA,r '(-ff) • < 50 > 



a, . in degrees 
dir 


Fig. 10. RMS errors of wind direction as a function of 
pressure and elevation angle. The values are 
maximum for the stated elevation angles and 
pressures. Results at 12.5° are based on data 
supplied by the Air Force Missile Test Center 
(1963). 



RMS errors for an elevation angle of 12.5° are given in Fig. 10 for 
the same scalar wind speeds that were used previously. Good agreement 
can be seen between these values and those cited for the AVE IIP data. 

The results presented above are not strictly valid when angle 
data are available at only 1-min intervals. These angle data are 
interpolated to 0.5-min intervals which constitutes a smoothing process 
on the input values. The degree of smoothing is dependent on the angle 
profile, but generally, the accuracy statements made earlier are still 
applicable . 

In conclusion, it may be stated that the rms errors cited above 
for the "raw" winds are in general agreement with those obtained from 
other sources and are highly dependent on balloon location and ambient 
wind speed. 

2. Accuracy of the Smoothed Winds . Figure 11 shows "raw" wind 
speed plotted as a function of time from release for Eglin Air Force 
Base, Florida at 1200 GMT on II May, 1974. In this example elevation 
angles ranged from 16° shortly after release to 30° near the top of the 
sounding. The figure indicates many small-scale variations in wind 
speed that become more apparent with increasing altitude. While some 
of these variations are indeed mesoscale meteorological phenomena, many 
are simply due to inaccuracies in the angle data which were truncated 
to the nearest 0.1° at 30-sec intervals. It was felt desirable to 
smooth these winds so that very small-scale features would be suppressed 

"Raw" wind speeds for Rapid City, South Dakota at 0300 GMT on 12 
May, 1974 are shown in Fig. 12 where the elevation angle ranges from 


Time From Release (min) 
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about 20° near the surface to 6,5° at about 40 min after release to 
12° near the top of the sounding. The wind speeds are observed to 
become very erratic at elevation angles below about 10° with most of 
these fluctuations undoubtedly being non-meteoro logical in nature. 

Similar results were obtained in the AVE X data when elevation angles 

. . ^ 

were less than 10°, Danielsen and Duquet (1966) found that errors 
in the elevation angles increase to 10 and 15 times the generally 
accepted rms error of 0.05° when the elevation angle is less than 10° 
They concluded that the increase in error was due to the inability of 
the tracking unit to discriminate between direct and ground -reflected 
or refracted signals. Some method of smoothing is clearly needed to 
suppress the large oscillations in wind speed associated with low 
elevation angles. 

Several smoothing techniques were investigated in order to find 
a procedure xvhich would perform the types of smoothing needed, be 
practical from a programming and cost viewpoint, and lend itself to 
error analysis procedures. Smoothing the angle databy least square 
polynomial fitting was considered, but results obtained when this 
procedure was used on AVE I .data were not satisfactory according to 
personal communication with NASA personnel. The application of 
smoothing functions to the elevation angles and finite differencing 
over larger time intervals were considered but not used. 

* Personal communication with Mr. W. W. Vaughan, Chief, Aerospace 
Environment Division, NASA Marshall Space Flight Center, Huntsville, 
Alabama, 35812. 
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It was finally decided to apply 5 -point binomial coefficients to 
the "raw" wind components that were, obtained at 30-sec intervals, and 
then to recompute scalar wind speed and direction using the smoothed 
components. The response function for the 5-point weighted average is 
shown in Fig. ,13. Changes in speed with a frequency higher than 1 cycle/ 
min are completely eliminated by the smoothing process while events of 
lower frequency are smoothed to a lesser extent. For example, 817, of 
the amplitude of changes in speed with a frequency of 1 cycle/5 min 
would be retained by the procedure, and 25% of the amplitude of changes 
in speed with a frequency of 1 cycle/2 min would be retained. Assuming 
an ascent rate of 1000 ft /min, this means that 81% of the measured 
shear in a 2500-ft layer of the atmosphere would be retained by the 
smoothing technique, and 25% of the measured shear in a lOOQ-ft layer 
would be retained. Angle data were available at 0.5-min (approximately 
500-ft) intervals. Based on an analysis of the wind data, it was 
decided not to begin the smoothing process until the balloon was 2.0 km 
above the ground. Figure 14 shows the results of the smoothing process 
which for purpose of illustration was not begun until 23 min after 
release. A comparison of Fig. 11 with Fig. 14 reveals the effect of 
smoothing. Tire accuracy of the smoothed winds is greater than, that of 
the original winds although an exact comparison of the accuracies cannot 
be made. 

3. Accuracy of Pressure Contact Winds . Smoothed wind components 
at 0.5~min intervals were interpolated on the basis of time from release 
to correspond to data given for the pressure contacts. New scalar wind 


Res 





Time From Release (min) 





speeds and directions were then computed. Near the surface, pressure 
contacts occur at approximately 0.5-min intervals so that no great 
amount of smoothing is involved in matching the winds with the contact 
data. At higher altitudes, the difference in time between successive 
pressure contacts increases so that at 200 mb the difference is about 
0.6 min, at 50 mb about 1.0 min, and at 25 mb about 1.5 min. At these 
higher levels, more and more smoothing is involved in the interpolation 
process. Figure 15 shows the results of the interpolation process on 
the smoothed winds in Fig. 14. The added smoothing due to the larger 
time interval between pressure contacts is evident near the top of the 
sounding. This is not undesirable because the elevation angle generally 
decreases with altitude causing poorer quality wind data. 

Figure 16 shows the results of smoothing and interpolation to 
pressure contacts of the wind speeds shown in Fig. 12 where the 
elevation angle becomes as low as 6.5°. A great improvement is evident 
although there are still fluctuations with periods of several minutes 
and amplitudes of several tens of m/sec. Figure 17 shows a time series 
of wind profiles for Topeka, Kansas, in which the minimum elevation 
angle for each profile decreased from 10.1° at the first time period 
to 6.8° at the last time period. Major features of the wind field can 
be tracked with time in spite of the decreasing elevation angle. 

Figure 18 reveals the same tracking capability in cases of higher 
elevation angles. Instead of resorting to even more powerful and 
complicated smoothing techniques, it was decided to mark wind values 
on the output with an asterisk if they were based on an elevation angle 



Fig. 15. Wind speeds on the pressure contacts at F,gl in Air 
Force Base, Florida at 1200 GMT on 1.1 May, 1974. 
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Fig. 16. Wind speeds on the pressure contacts at Rapid City, 
South Dakota at 0300 GMT on 12 May, 1974. 




Time From 



Wind Speed (m/sec) 

pig. 17, Time series of wind speed at Topeka, Kansas from 1200 GMT on 11 May, 19/4 
to 0000 GMT on 12 May, 1974. 





Time From Release (mm; 


20.4° 27.5° 22.3° 18.0° Min. Elevation Angle 

11/15 GMT 11/18 GMT 11/21 GMT 12/00 GMT Time 



0 5 10 0 5 10 

Wind Speed (m/sec) 

Fig. 18. Time series of wind speed at Charleston, South Carolina from 
1200 GMT on 11 May, 1974 to 0000 GMT on 12 May, 1974. 
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that was less than 9°. This action seems justified based on Figs. 17 
and 18 and the. accuracy of data presented earlier for winds based on 
low elevation angles. Users can then determine their own use of the 
particular piece of information. It should be noted that there are 
relatively few soundings which contain elevation angles as low as those 
described in Figs. 12, 16, and 17, 

The accuracy of the AVE IIP (average) wind data on pressure contacts 
is greater than that stated for the "raw" winds because of the added 
smoothing and interpolation of this stage. In addition, errors cited 
for the "raw" winds were maxima for the stated elevation angles and 
pressure surfaces. 

4. Accuracy of the Winds at 25-mb Intervals . Linear interpolation 
was used to obtain winds at 25-mb intervals. The winds being interpolated 
had previously been smoothed and interpolated to pressure contacts as 
described earlier. The accuracy of the winds is similar to the values 
cited p /iously. 

C. Comparison of AVE IIP Results with that Obtained by the National 
Weather Service 

It is desirable that the results of the AVE IIP reduction process 
be compatible with those obtained by the National Weather Service so 
that, if required, data from both sources may be used interchangeably. 

This section will not seek to demonstrate the superiority of either 
reduction procedure although the pressure contact results that are 
derived from AVE IIP data contain far more information than that 
received from the NWS. Figures 19-21 contain sounding data processed 
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Fig. 20, AVE IIP and National Weather Service soundings for Flint, 
Michigan at 1200 GMT on 11 May, 1974. The minimum antenna 
elevation angle was 12.0°. 
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Fig. 21. AVE IIP and National Weather Service soundings for St. Cloud, 
Minnesota at 0000 GMT on 12 May, 1974. The minimum antenna 
elevation angle was 10.3°.' 
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by the AVE II procedure and data that is transmitted nationally over 
teletypewriter circuits by the NWS. Figures 22-35 show constant-pressure 
charts obtained from either AVE IIP data or NWS data for the 850- , 700- , 
500-', 400-, 300-, 200-, and 100-mb levels for 1200 GMT on 12 May, 1974. 

The two sets of figures will be the basis of comparison between the two 
methods of data reduction. 

1. Temperature Comparisons . The temperature computation scheme 
used in AVE IIP is the same as that currently used by the NWS and 
previously described in Section II-C of this report. A comparison of 
Figs. 19-35 reveals a very close agreement between the two sources of 
data. The similarity is especially evident in the three soundings. 

2. Humidity Comparisons . The procedure used by the NWS for 
computing humidity is somewhat different from that used in AVE IIP 
(Section II-C). The NWS system uses a baseline ordinate corresponding 
to a relative humidity of 50%, after the baseline data has been locked 
in on their evaluator. Their reduction equations make use of constants 
in all computations with the constants being derived from the baseline 
check. Separate equations, depending on whether high or low ordinate 
values are encountered, are used in the evaluation process. Different 
procedures are used by the NWS to evaluate the dew point temperature 

as well. The NWS procedure sets humidity at 20% if the computed value 
is less than 20% while the AVE IIP procedure outputs the computed value 
down to a value of 5% where it is then defined to be missing. Dew point 
depression is set at 30C at relative humidities lower than 20% in the 
NWS procedure while the AVE IIP scheme computes values until the relative 
humidity is below 5%. 





















Fig. 28. 400 -mb chart for 1200 GMT on 12 May, 1974 using AVE IIP 
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Fig. 32. 2 00 -mb chart for 1200 GMT on 12 May, 1974 using AVE IIP 
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These differences in computational procedure appear as slight 
differences in dew point temperature in Figs. 19-21; the difference 
in data detail is quite noticeable. The constant pressure charts shown 
in Figs. 22-35 show slight differences in the results of the two 
procedures as well. No large systematic differences between results 
of the two procedures are indicated although NWS humidity values appear 
to be slightly lower than AVE IIP values in some cases, especially at 
low values of relative humidity. No significant problems should be 
encountered, however, in using NWS data around the borders of the AVE 
II pilot experiment area. 

3. Height Comparisons . The same procedures for computing height 
are used in both the NWS and AVE IIP reduction schemes. Many more data 
points are used in the AVE IIP procedure so that these results should 
be superior to those obtained by the NWS. Figures 22-35 reveal only 
slight differences between heights, and there does not appear to be 
any systematic difference between results of the two procedures. 

4. Wind Direction and Speed . The wind reduction procedure of the 
NWS differs from that used in AVE IIP which is described in Section II-C 
of this report. The NWS receives angle data at 1-min intervals instead 
of 30-sec intervals as obtained in the AVE II pilot experiment. The NWS 
scheme involves finite differencing over 2 -min intervals below 14 km 
and over 4-min intervals above 14 km. If the elevation angle is less 
than 12°, it is smoothed by use of a 3-min average before the differencing 
process begins; the reduction process terminates at elevation angles 
below 6°. Figures 19-22 reveal no large differences between results of 
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the two procedures, even in areas of high wind speed such as in Fig. 19 
for Peoria, Illinois. 


In conclusion, the results of this section indicate that matching 
AVE IIP and NWS values of temperature, moisture, height, and wind at 
the borders of the AVE II data should pose no problems. 
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APPENDIX A 

General Features of the Reduction Programs 

i 

The reduction programs used at Texas ASM University were written 
using WATFIV, a FORTRAN IV compatible language developed at the 
University of Waterloo, Ontario, Canada. WATFIV is almost completely 
compatible with FORTRAN G and H, but the CHARACTER * declaration state- 
ment used in some of the programs cannot be used in FORTRAN G or H. 

Other restrictions between the various compilers may be found in texts 
on computer programming. 

The data reduction process at Texas A&M University was performed 
using an IBM 360/65 computer with OS/MVT (Multiprogramming with a 
Variable Number of Tasks) and HASP (Houston Automatic Spooling Program). 
Use of these programs on other machines may require considerable program 
modification. 

The comments made in the following appendices are meant to 
supplement the program descriptions that have previously been given and 
are mainly concerned with required input data and core allocation for 


the programs. 
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APPENDIX B 

Use of the Card to Tape Program 

Card decks described in Tables 2-4 are transferred to nine track 
tapes that are non-labeled. The DCB (Data Control Block) parameters 
are DCB » (RECFM = FB, LRECL = 80, BLKSIZE = 1600, DEN =3). An 
object code of 5,608 bytes and an array area of 40,160 bytes are 
required. The data required for the program are one card specifying 
the number of soundings to be processed followed by the stated number 
of complete sounding decks. No printed output is produced by the 


execution of this program. 




orvo noo oon ooo non non r> o n o ooo non o o o 
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C TRANSFER THE CAROS TO TAPE 

C SUMMFR 1 97 A — HENRY FUELBERG 

C CHARACTER * IS AN IMPORTANT FEATURE OF WATFIV 

C 

CHARACTER*! CMB.OJ , IOC 80), ONE/* l •/ 

char act er* 8 q card! 5oo>,cardi 

EQUIVALENCE ICI, CAROI) 

REWINOl 

READ THE NUMBER OF SOUNDINGS TO TRANSFER 

READ , NSQUND 
00 2 0 M= !, NSOUNO 
1=1 

READ THE CARD 

4 RFAD (5 , 100 ) CARDT 
100 FORM AT (A 80) 

IF ( I .NE. 1) GO TO 7 

COMPUTE THE NEW 'LEADER CARD* -USING BASELINE INFORMATION 

DO 6 L = 50, 69 
LL=L-49 

6 ID ( L L) =C H L ) 

7 CARD (I )=CARDI 

IF l I .EQ. 1) GO TO 2 • 

A *1* IN COLUMN 72 INDICATES THE END OF THE ORDINATE DECK 
A *1' IN COLUMN 72 OF THE BASELINE CARD IS NOT CONSIDERED 

• IF 1C II 72 ). EQ.ONE ) GO TO 3 

COUNT THE NUMBER OF CARDS IN THE ORDINATE DECK 

2 1 = 1*1 
GO TO 4 

ADD THE NUMBER OF CARDS IN THE ANGLE DECK 

3 1=1 + 1 

READ THE CARD 

RE AO ( 5 , 1 00 ) CAROI 
CARD (I )=CARDl 

A *1* IN COLUMN 72 INDICATES THE END OF THE ANGLE DECK 

IFIC It 72). EQ.ONE) GO TO 10 
GO TO. 3 

WRITE THE 'LEADER CARD* ON TAPE 

10 WRITE! 1, 5) ( ID(LL) ,LL=1,20),I 

5 FORM AT (20A 1,14) 

TRANSFER THE REMAINING CARDS TO TAPE 
WRITE! 1 , 15 ) (CARD! J ) » J = 1,1 ) 


15 FORM AT ( A80 ) 
20 C0N7 INUE 
REWINOl 
STOP 
END 
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APPENDIX C 

Use of the Tape to Print Program 

No input data cards are required since the tape is read and printed 
until an END OF FILE mark is encountered. The same tape parameters that 
were used in creating the tape are used to read from it. The object 
code is 4,864 bytes while the array area is 80 bytes. The printed 
output resembles the original data decks. 


C THIS PROGRAM PRINTS THE RAW DATA TAPE 

C SUMMER 1974— HENRY FUELBERG 

C 

DIMENSION NCARD120) 

REWIND! 

ID READ <1 . 11. END=13> (NCARD(I),I=l,20) 

U FORMAT (20A4) 

WR ITE! 6. 12} (NCAROi I > . 1*1, 20 ) 

12 FORM AT ( IX, 20A4) 

GO TO 10 

13 REWIND! 

STOP 

END 
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APPENDIX D 

Use of the Master Reduction Program 

Angle, ordinate, and baseline data are read from magnetic tape 
with the DCB parameters described in Appendix B. The input format is 
described in Tables 2-5. The output data is placed on a non- labeled 
tape with DCB ® (RECFM = VS, DEN =3); no length of record or block size is 
specified. Data read from cards include the number of stations in the 
station roster, in this case 54, followed by a card describing each of 
the stations. The last card read specifies the number of consecutive 
soundings to be processed. The object code of the program is 29,840 
bytes with an array area of 39,136 bytes. Printed output from the 
program is described in Section II-C and Fig. 5; the identical output 
is transferred to magnetic tape and is described in Table 7. The sounding 
identification data is followed by a series of 18 arrays in the order 
given in Table 6 that correspond to the contact data. Each array has 


dimension of 230. 


o o o o n o o o o o o o o o o o o cv n o o o orv o o o o o r> o 


C THIS PROGRAM REDUCES RAH RAWINSINDE DATA TO A FINISHED PRODUCT. 

C SUMMER 1 07 A HENRY FUELBERG 

C SONOES ARE WEATHER SERVICE TYPE WITH CARBON HUMIDITY ELEMENTS 

C Cl PENSION NECFSSARY ARRAYS 

C 

DOUBLE PRECISION DEXP 

D I MENS ION R I 2 BO) ,ATH(230> , XX S( 230) ,ZZS(230) 

DIMENSION WDD( 230) ,WSS (230) , CD (5 ) 

DIMENSION T MORI 2 30), P I 2 30 ) , TQRD { 2 30 ) , HGRDt 230) , TMNG(270), 

1 THETAt 270) , AZ( 27 0) , XSC 230) ,ZS(23 0) .HMC230) , TC < 2 30 ) , WD < 230 ) , 
2VS1230 ) ,C( 40), KC(20) ,TVt 230) , TD( 230) ,W( 2 30) , PTK( 230 ) .HINT! 230) , 

3 STHT ( 6 0 1 , S TLA [ 60 ) , S TLG ( 60 ) . EP OT (2 30 ) , PH I ( 23 0 ) ,CTC(230> 

DIMENSION YS(230), WWEt 230) ,KSN(230 ) , WE (230) , WN(230) 

INTEGER ST I 0(6 0) , [STOP , A STOP ,1 ASMt 27 0) , I E P ( 2 70 ) , I ORI N( 230) 

INTEGER 8LANK.AST 

INTEGER IFP 1230) ,NAME( 7) ,NMS 170, 7) ,NZ<6) 

DIMENSION U ( 230 ) ,Vt 230) 

RFSISTITORD) = DEXP ( 16 . 00 8299L-0.9966256«AL0G( 2*T0R0 >) - 48000. 
XVIRT1 XTMP , XE,XP1=XTM? /( I . 0- l 0.3 79* < XE/XP) )) 

C 

C DEFINE HUMIDITY CONSTANTS TO BE USED 

C 

DATA C/122 0. 7,-19. 55,8.3 12,-. 12738 ,'1249. 8,-1 9.307, 11 .7 32 . 1 7732 , 
L -1687.8, 26.819,-6. 1439 ,.096 15, -201 5. 2, 30. 6 85, -12. 115*. 17848 

2 34 1V2 » — 5. 8167, -1.9779, .024 92, 7 74.05,-1 1. 77 3, 1. 10 14 . 0 1397, 

3 22. 327, -.4 0763,. 31303,-. 00344,-112.37, l. 67 12, -.04262 , 

4 -.00017,-57.854, . 5 39 14 I 5883 ,.00234, 8* 0627 97 356, 

.5 . 13345, . 00032/ 

CATA A S T , 3 LANK/ IHM , 1H / 

RD=6.857E-2 


DEFINITION CF IMPORTANT VARIABLES 

ASTOP THE LAST CARD IN THE ANGLE DECK IS INDICATED BY A •l* IN 
COLUMN 72 

ATH IS AZIMUTH OF SONDE MEASURED CLOCKWISE FROM NORTH 

AZ IS ELEVATION ANGLE 

CP IS SPECIFIC HEAT OF MOIST AIR 

CTC IS PRESSURE CONTACT NUMBER 

OT IS THE INTERPOLATION WEIGHTING FACTOR 

E IS VAPOR PRESSURE 

EL IS LATENT HEAT OF EVAPORATION 

EPOT IS EQUIVALENT POTENTIAL TEMPERATURE 

HINT IS HEIGHT AT THE ANGLE TIMES, 30 SEC APART 

HM IS RELATIVE HUMIDITY 

HORD ’ IS HUMIDITY OR 0! NAT E 

TlMIN IF EQUAL TO 0 ANGLES ARE EVERY 30 SEC 

IF EQUAL TO l ANGLES ARE EVERY 1 MIN 

I ASM A *1' IN COLUMN 73 INDICATES THAT MANUAL SMOOTHING WAS 
DONE DN THE DATA CONTAINED ON THAT CARO 
ID! IS THE DAY DF THE SOUNDING 

ID2 IS THE MONTH OF THE SOUNDING 

103 IS THE YEAR OF THE SOUNDING 

!D4 IS THE TIME OF THE SGUNOING 

IDS IS THF STATION NUMBER 

I E.P A >1' IN COLUMN 74 INDICATES THAT THE ELEVATION ANGLE IS 

LESS THAN 9 OG. THESE ARE ALSO CHECKED BY THE PROGRAM 
IFP IS- USED TO MARK WINDS THAT WERE COMPUTED W I TH .ELEV AT ION 

ANGLES LESS THAN 9 DG 

IORIN A ' l* INDICATES THAT VALUES CN THIS CARD HAVE BEEN 
INTERPOLATED 

ISTOP A '1' IN COLUMN 72 TO INDICATE THE END OF THE ORDINATE 


o o anr% o >*» v> r> o n o n no n o o n o o o o o o o o n r> o o o o r* o o o o o o r* o o 
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DATA 

P IS ATMOSPHERIC PRESSURE 

PTK IS POTENTIAL TEMPERATURE 

R IS RANGE OF SONOE FRO* RELEASE 

RESIST IS A FUNCTION USED TO COMPUTE TEMPERATURE 

RHO IS HUMIDITY AT ORDINATE 46 ANO TEMP«-40C 

PWSFC IS SURFACE RELATIVE HUMIDITY 

SC IS THE SONDE'S DISTANCE OVER A CURVED EARTH 

TC IS TEMPERATURE 

TO IS DEW POINT TEMP 

THETA IS ELEVATION ANGLE 

TM I IS MEAN VIRTUAL TEMPERATURE 

TMNG IS TIME OF THE ANGLE OBSERVATION 

T M OR IS TIMF OF THE CONTACT 

in IS THE BASELINE TEMPERATURE AT 37.6 OROINATES 

TO HAS AN IMPLIED NEGATIVE VALUE 
TPRD IS TEMPERATURE ORDINATE 

T SC. IS THE APPROXIMATE TEMP AT THE LCL 

TSFC IS SURFACE TEMPERATURE 

TV IS VIRTUAL TEMPERATURE ' 

U IS THE EAST-WEST WIND COMPONENT AFTER INTERPOLATION 

V IS THE NORTH-SOUTH WIND COMPONENT AFTER INTERPOLATION 

W IS MIXING RATIO IN GM/KG 

WO IS WIND DIRECTION 

won IS WIND DI RECT ION 

WO I IS SURFACE WIND DIRECTION 

WF ts THE EAST-WEST WINO COMPONENT AFTER SMOOTHING 

WN IS THE NORTH-SOUTH WIND COMPONENT AFTER SMOOTHING 

VS! IS SURFACE WIND SPEEO 

VS IS SCALAR. WIND SPEED 

VSS IS SCALAR WIND SPEED 

WWE IS THF EAST-WEST WIND COMPONENT BEFORE SMOOTHING 

WSN IS THE NORTH-SOUTH WIND COMPONENT BEFORE SMOOTHING 

X IS MIXING RATIO IN GH/GM 

XS IS THE X LOCATION COMPONENT 

XV I R T IS a FUNCTION USED TO COMPUTE THE VIRTUAL TEMPERATURE 
XxS ts THE X LOCATION OF THE SONDE ON A PRESSURE CONTACT 

YS IS BALLOON HEIGHT 

l S IS THF Y LOCATION COMOONENT 

ZZS IS THE Y LOCATION OF THE SONDE ON A PRESSURE CONTACT 

REWIND1 

REWIN03 

DEFINE ALL ARRAYS TO BE PRINTED LATER 

DO 3 45 1=1.230 
TMOR (I ) = 0.0 
CTCm=0.0 
YS< I >=0.0 
PI I > =0.0 
TCI I 1=0.0 
lORIHI n=o 
roil 1=0.0 

VDDI 11=0.0 

wssi n=o.o 
ifpii)=o 
ut 11=0.0 ■ 

VII»*0.0 
PT KI 11*0.0 . 

EPOT II I *0. 0 


i. 

i 
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m n =o.o 

HMl ! >=Q.O 
R( II =0.0 
ATH( I ) =0.0 
345 CONTINUE 

r. 

C READ TUESTATION ROSTER CONTAINING WHO NUMBER, STATION ELEVATION, 

C STATION LATITUDE, ANO STATION NAME 

C NR IS THE NUMBER OF STATIONS IN THE ROSTER 

C 

READ, NR 
DO 111 J=l , NR 

READ (5, SOU -STIOCJI.STHTI j ) , STLA { J ) , (NMS ( J,K > , K= 1, 7» 

501 FORMAT { 1 5 , 1 X , F4 .0 , IX, F A. 2 , 2 AX , 7 A4 ) 

111 CCNT INUE 

RE AD.NSOUND 
C 
C 

C N SOUND IS THE NUMBER OF SOUNDINGS TO COMPUTE 

C EFGIN THE GRAND LOOP 

C 
C 

00 9 99 I S= UNSOUND 
C 
C 

C RE AO THE BASELINE DATA 

C THE FIRST BASELINE ’CARO' NEED NOT BE READ UNLESS ONE INTENDS TO 

C SEARCH THE TAPE. THE TOTAL NUMBER OF CARDS IN THE SOUNDING DECK 

C IS CONTAINED ON THIS. CARD 

C 
C 

RFAD(l.lO) <NZ( I), 1=1,6) 

10 FORMAT (0^40 
C 

C READ THE SECOND 'CARO 1 WHICH CONTAINS BASELINE CALIBRATIONS AND 

C SURFACE MEASUREMENTS 

C 

RE AD (1,274 , ERR =9 36 61 C TC (1) , PC 1) , T SF C, RH SFC , TO, RHO , WS I , WDI , 

: l ID2, 101, ID 3, 104, 1.05, UMIN 

2 74 FORM AT ( 7X, F4.1 , IX, F6.1 , 1X.F4.1 , 1 X, F4.1 , 1X,F4.1 ,1X,F4.1, IX, F4.1, 

1 1X,F3.0,2X, 12, IX, I 2, IX ,12, IX, l 4, IX, I5.3X, II) 

TM0R(1)=0.0 

C 

C THE BASELINE TEMP (TO) CORRESPONDS TO AN ORDINATE (TCAU OF 37.6 

C 

TO =- 1. 0*TO 
TC AL=3 7. 6 
IOR I N< 1 )=0 
C 
C 

C SET UP THE BASELINE REFERENCE 

C WEATHER SERVICE TEMPERATURE ELEMENT 

C 
C 

RT B= RES 1ST ( TCAL ) 

RK l=TO + 273 .15 

RM 1=5. 30.1B981X 1 .0/ 303 . 0-1 . 0/R K1 J 

RM 2= (-2.47991E-3+SQRTC 2. 4799 IE-3 ** 2-4. *5 . 899 86E-5*RM 1 > )/ 

1(2.* 5. 8993 6E-5) 

RM3= 14000. /RTB*EXP(RM2 » 

C 



107 


C SFT UP CONSTANTS TO BE USED IN HUMIDITY CALCULATIONS 
C CARBON HLM1DITY ELEMENT 

C 

KX=0 

DO 103 j =1,40, 2 
KX = KX ♦ l 

103 FC(KX) = CIJI * 1C I J*l l*RHOI 
C 

c 

C READ the ORDINATE OATA 

C 

C 

DO 233 10=2,230 

PE AO (I , 234 , ERR =0366 I C TC I 101 ,P ( I 01 ,TOROI 101 ,HOROI|OI ,TMOR| 101* 

1 ISTOP, ! OP INI 101 
IH 1 STOP. NE. 01 GO TO 235 
233 CONTINUE 
235 KL M = 10*1 

2 34 FORMAT I F 3. 0 »1X »F 4. 0 . IX ,F 4. 1 , IX ,F 4. I , IX , F5 . 1 , 47 X, ■ I , I 1 1 
C 

C KLM IS THE NUMBER OF CONTACTS 

C CHECK ORDINATFS FOR MISSING TIME VALUES 
C USE LINEAR INTERPOLATI CN TO FILL IN GAPS 

C PRESSURE IS ALWAYS GIVEN F VEN |F TIME IS NOT GIVEN 

C 

LT IN =0 
IT J= 0 

DO HU 1= l.KLH . 

IFILTIN.EO.OI I T J= 0 

IF ( T MOP.l II . NT • 999, 9 ) GO 'i'O l III 

IF 1 TMORC II . El!. 994. 4 ) T 1= THOR 1 1 -l I 

IF I LT IN.NF .01 GO TO l II 5 

DO I II 2 J*l»KtM 

IT J= |T jh 

IFtTMORI Jl .NE.999.91 GO TO 1113 

1112 CONTINUE 

1113 T 2 = T MO R I J ) 

DT = T 2-T l 
8DT=CT/1TJ 

1115 T M QR It I = TMOR I l-l >» AOT 
TOPOU 1 = 99.9 
HOROII 1=99.9 

r 

C INDICATE THAT INTERPOLATION WAS DONE 

C 

ICR INI I »=l 
LT IN » I TJ-t 
till CONTINUE 

IF I T1MIN.NE.0I GO TO 3333 
C 

c 

C READ THE ANGLE DATA IF HALF MINUTE VALUES 

C I r TMNG (ll.NE.O, INTERPOLATE THE TIMES BUT FILL IN THE MISSING 

C ANGLES WITH 9*S 

C 

C L -,. - . ■ ■ 

CO 237 t A= 1 , 270 

READ! I ,2 33 ,FRR =93661 T MNG 1 1 » I , THEY A 1 I A I , A7.I I A) • AST OP « I ASH 1 1 A! , 
1 IEPI tAI . ' . - 

IFITMNGI IA I .EQ.999.9) GO TO 3777 
INC= (TMNG1 n/o.Sno. I 






I 
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T M NG ( I NGt L ) = TMNG I I A ) 

THETAl ING+- 1 ) =THETA( IA) 

AZ< INCU)=AZ( I A ) 

IASM UNG+l ) = I ASM t I A > 

1 F P ( INCH) = I 6P < IA) 

time=q.o 

IF( I NG.eO.O ) GO TO 3454 
DO 237a IZ = 1 , ING 
T M NG ( I Z)=T IME 
THETAl I2)=99.9 
IFP(IZ)=0 
AZ ( I Z ) =999 . 9 
I ASM ( I Z )=0 
T IME=T I ME«-0 • 5 
2378 CONTINUE 
GO TO 3454 

3777 IFl ASTOP.NE.O) GO TO 239 

237 CONTINUE 
3454 I NG = I NG+2 

C 

C READ THE ANGLES IN THE NORMAL FASHION NOW 

C 

CO 2377 I A = I NG» 2 70 

READ ( 1 ,238 , ERR=9366 ) T MNG ( I A ) , TH ET A II A ) , AZ (I A) ,A ST OP , I ASM ( I A) , 
1 IE P ( l A ) 

IF (ASTOP.NE.O) GO TO 239 
2377 CONTINUE 

238 FORM AT (F5. I , IX, F4. I , IX ,F 5. 1 , 55 X, U , 1 1, 1 1 ) 

239 K A T = I A — I 
IF (KAT.GT. 230) KAT=228 
GO TO 4444 

READ THE ANGLE DATA IF WHOLE MINUTE VALUES 
3333 00 5555 IA =1,270, 2 

READ! 1,2388, ErR = 9366) TMNG ( l A) , THE TA ( I A ) , AZ (I A ) , A STOP , I AS M (I A ) 
llEP(IA) 

TMNG(IA«-1) = TMNG< IAJ4-0.5 
THET A I I A + l ) =99.9 
A Z ( 1 A+l>=999.9 
I ASM ( t A t l) ”0 
16 P I IA + I)=0 

IF( ASTOP.NE.O) GO TO 2399 
5555 CONTINUE 

2388 FORMAT (F 5. L,1X,F4. L, IX ,F 5. 1 , 55 X, U , 1 1, 1 1 ) 

2399 KAT= IA-1 

IFIKAT .GT.230) KAT=228 
C 

C INTERPOLATE WHOLE MINUTE ANGLES TO HALF MINUTE VALUES 

C 

DO 6666 K= 2 » KA T * 2 

IF < THETA (K + D.EO. 99.9) GO TO 6667 
IF(THETA(K-n.EQ.99.9> GO TO 6667 
IF(AZ(K+1).E0. 999.9) GO TO 6667 
IF < A Z ( K-l) .ED. 999. 9) GO TO 6667 

THFTAt K)=THETA (K-L) +0. 5* C THE TA ( KH ) -THET A( K- 1) ) 

CAZ=AZ(Ktl )-AZ(K-l) • 

AZ (K )= AZ(K- HtO. 5*DAZ 
I F ( A3S ( OAZ ).GT .180.) A Z ( K ) = A Z( K ) *1 80 . 

I F ( A Z ( K ) . G E • 3 6 o . ) AZ (K ) = AZ< K 1-360. 

■ .GO TO 6666 
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6667 THET MK!*99.9 
A7 IK 1*999. 9 
6666 CONTINUE 
C 

C CHECK THE ROSTER FOR THE PARTICULAR STATION TO ASSIGN 

C ELEVATION A'NO LATITUOE 

C 

9996 00 HZ J=l,NR 
JS = J 

IFIID5.EQ. STIPI Jll GO TO 113 

112 CONTINUE 
WRITE! 6,96 6 1 

966 FORM at I IMl , 39HERR0R — S T AT ION 10 IS NOT IN ROSTER I 
GO TO 999 

113 HT = S TUT I JS I 
00 26 1 = 1, 7 

NAME ( I ) = NM S IJS , I I 
26 CONTINUE 

1F( I05.EQ. 259) 105=260 

AZIMUTH ANGLES AT NSSL STATIONS ARE 180 DEG OFF 

IF II 05. LT. 220021 GO TO 1213 

IFIID5.GT. 22005) GO TO 1213 

Pn 1 212 NS S L = I , K AT 

IFIAZINSSL I .EQ. 999.91 GO TO 1212 

AZ ( N SSL I =AZ I NSSL !♦ I BO. 0 

IFIAZINSSL I .GE .360. 0.1 AZ I NSSL) = A Z( NS SL I- 360. 0 

1212 CONTINUE 

1213 CONTINUE 

CONVERT NUMERICAL MONTH TO ALPHABETICAL 

CALL MOPPT ( 102, IX, I V, I Zl 
TCI l 1= T SFC 


COMPUTE TEMPERATURE AT THE PRESSURE CONTACTS 
WEATHER SERVICE TEMPERATURE EQUATIONS 


00 999 K=2,KLM 

IFITOROfM .EQ.99.9) GO TO 993 
RE = PES IS T C TORO IK I ) 

R9= A LOG IRM3*REZ19000.0I 

P.K=I./( l •/ 3 03. *-9.67 79E-9*R9* 1. 1 12TBE -5*R4**2I 
TP-RK-273. IS 
TC I K )= TP 
GO TO 999 
TCIK 1 = 99.9 
CONTINUE 

CHECK FOR MISSING TEMPERATURES 
USE LINEAR I NTERP0LAT1 CN TO FILL IN GAPS 

00 991 1 1=1, KLM 
IF IT Cl II.NE.99.9I GO TO 9.911 
00 9612 J* 1 ,KLM 
KOJ= J 

IF ( TORDt Jl . NE.99.9 1 GO TO 9913 
9912 CONT INUE 


non opo oono oooo ■ .00=00 .000 r> o o 
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*<►13 OTKT=(TMOU ! I-TMORI 1- 1 ) ) / l TMOR l KG J ) -THOR ( I- 1 ) ) 

TC (I )=TCtI-l>+DTKT*lTC (KOJ)-TCU-l)) 

IORINI T > = 1 
4411 CONTINUE 
C 
C 

C COMPUTE OTHER THERMODYNAMIC VARIABLES AT THE CONTACTS 

C 

C 

00 A A K=l, KLM 
TP=TC( K t 
TK=TC( KJ+273.16 
IF (K .EQ. 1) GO TO 787 

IF HUMIDITY ORDINATE EQUALS 99.9, PRINT 9«S FOR MOISTURE VARIABLES 

IF(HORDtK) .EQ.99.9) GO TO 523 
IF (HGRDIKl i 523, 523, 52 

X AS USED HERE IS NOT MIXING RATIO 

52 X = IHORDl K )-A6. )/ Al. 

HP = 0. 

DO 521 KK= 1,17,2 

52 L HP = < HP+HC (K'<) + HCtKK+ll+TP) * X 
HP = HP + HCI19I + HC( 20 ) * TP 

MOISTURE VARIABLES BASED ON A RELATIVE HUMIDITY LESS THAN 5 PCT 
AP.E PRINTED AS 9*S 

522 IP I HP-5. ) 523, 532, 532 

523 HP=999.9 

TV t K ) = TK 
TD(KI - 99.9 
E = 0.0 

PT K t K ) = TK* ( t 1000./ (P IK)) )**.286> 

EPOT (K )= 99 9 . 9 
W! K) =99.9 
HM(K)=HP 
GO TO 537 
787 HP =R HSFC 

THIS SECTION COMPUTES THERMODYNAMIC VARIABLES WHERE HUMIDITY IS 
GREATER THAN 5 PCT 

532 HM (K )=HP 

IF I HP- 100. ) 53A, 53A , 53 3 

PRINT HUMIDITY EVEN IF GREATER THAN 100 
OTHER VALUES BASED ON A MAX VALUE OF 100 

533 FP = 100. 

53.4 E = HP * .0611 * 10. **I ( 7.5*TP ) / (237.3+TPJ ) 

X USED HERE IS MIXING RATIO 

X= ( . £2 3* E) / t P I K ) — E J • 

CP=.2A0+(1 .0+0. 8 4* X) 

THIS WILL GIVE DEW POINT 
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ill 


FLOG = AL0G10IF) 

TD(K| * U Z37.3*ELUG)~186.527) /I8.206-ELOG) 

0= Tf IK l-TI) ( K1 

TSC= TU IK >- I 0.212*0.001 5U'*TD IK 1-0. 000*36*TCI K> 1*0 
I SA- T5C*27 3.16 
FL=597.3~. 506* TSC 

PTKI KMTK* ( ( 1000. /I P IK *-E> t**(RD/CP) I 
E*»OHK)«PTMK>*fcXPl IEL*X)/ICP*T$A|| 

WIKI = X*‘l 00 0 . 

TVIKI - XV IRTITK.E.PIKI) 

5 37 CONTINUE 

53 IMK-1 |61, 61,62 

C 

c 

C COMPUTE HEIGHT AT EACH CONTACT - 

C 

c 

61 S YS= 0. 0 
GO TO 63 

62 T M -I = I TV( K)>TV( K-ll I 72. 

DYS= 29 * 290 8 5*TM I *ALOG( PI K-ll /PlKll 
SY 5= SY S*DY S 

DON'T COMPUTE WINDS IF LESS THAN 10 ANGLE REPORTS 
INDICATE THIS HY DEFINING VALUES TO BE PRINTED WITH 9* S 

63 VS I K I = SYS 

IFIKAT.LT. 10) WSSt K I =9 9.9 
TF IKAT.LT* 10) WDDI K ) *999.9 
IFIKAT.LT. 10) IFP(K)*0 
IFIKAT.I.T. 10) I A S M I K | » 0 
lFiKAT.LT. 10) OIK) =99. 9 
1FIKAT.LT. 101 V I K I = 99. 9 
I r IK AT .LT. TO) R I K ) * 999 .9 
ir (KAt.LT. 10) ATHIKM999. 

A A CONT INUF 

IFIKAT.LT. 10) GO TO 921 


INTERPOLATE HEIGHTS TO CORRESPOND TO ANGLE T IMES 130 SECI 
LINEAR INTERPOLATION BASEO ON TIME TS USEO 


I J T = 1 

HINT 1 1 1 = 0. 0 
DO 701 1=2, KAT 

IFITMNGI I) . GT. TMOR I KLM1 I GO TO 77T 
IT- 1 

DO 702 1 J = I JT , 300 

IF ( T MOR ( IJ*1).GE.TMNG( 1)1 GO TO 703 

702 CONTINUE 

703 Tl= TMOR l IJ*1) 

T2=TM0R( IJ) 

DT=( TMNGlI >-Tl)/lTl-T2 ) 

HINT (I )= I YSt I J* II- YSI !J))*OT*YS(l J+l) 
I J T = I J 
701 CONTINUE 
777 KT*I T— | ■ 


COMPUTE WINOS AT 30 SEC INTERVALS. 1 NIN OVERLAP 
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c 

c 

00 750 K = 1 , KT 
C 

C CHANGE ANGL'ES TO RADIANS 

C I F AN ANGIE TO 06 USED IS HISSING (9'S), WIND VARIABLES AND SONOE 

C LOCATION WILL BE PRINTED AS 9'S 

C 

IF(THETA(K) .£0.99.9) GO TO 798 
IF(AZ( K) .EQ.999.9) GO TO 798 
THETA! K> =TH6TA ( K I / 5 7. 2 9578 
AZ(K)=AZ(K 1/57.29578 

As CO'SI THETA(K) ) / < l . ♦( H INT < K ) /6 37 12 29. ) ) 

SC = 6 37 12 29 .“(A RC OS ( Al-THETA! K) 1 
GO TO 799 

798 XS(K 1=99.9 
ZS (K 1=99.9 
IFCK-31 67,68,69 

799 XS(K>-SC*S INtAZIKl 1 
ZS(K)=SC*COS(AZ(K) 1 
IF (K-3 167, 68,69 

68 DT = 3 0. 

C 

C WIND AT THE FIRST LEVEL (21 ABOVE GROUND IS A. 30-SEC FORWARD 

C DIFFERENCE USING LEVELS 3 AND 2 

C 

XS(1 ) = X S ( 2 1 
Z S ( l ) = Z S ( 2 1 
GO TO ‘71 
69 DT=60. 0 

71 I F ( X S ( K ) • E 0 . 99 • 9 ) GO TO 751 

IF!XS!K-2) .£0.99.9) GO TO 751 
IF(ZSl K) .EO.99.9) GO TO 751. 

IF tZ SC K— 2) .EO.99.9) GO TO 751 

C 

C COMPUTE THE WIND COMPONENTS 

C 

WWE! K-l )={ XS!K)-XS( K-2 ) l/DT 

WSN( K-L )=( ZS (K)-ZS(K-2 ) ) /DT 

IF ( A BS ( WSN ( K-l ) ) .GT.l.OE— 51 GO TO 5151 

TT=WSN(K-I ) 

WSN ( K— 1 1 = l . 0 E— 5 

IF t TT. LT.O.O) WSNt K-l) = - 1. 0* WSN ( K- 1) 

5151 WS( K-l != SORT! ( KWE! K-l) 1**2 «• 1 W SN IK-1 1 1 ** 2 1 
8=WWF( K-l) /WSN(K-l) 

A= AT AN (A BS (8)1*57. 29 57 8 
I F ( B ) 671,671,672 

671 IF (WWE (K-l) 1 673,673,679 

673 A=360.-A 

GO TO 676 
679 A= 1 8 0. — A 
GO TO 676 

672 IF (WWE (K-l) ) 67e,678,676 
678 A= A + 18 0. 

676 CONTINUE 
C 

C MODIFY WIND DIRECTION TO THE GENERAL MET CONVENTION 

C 

A = A * 180. 

IFfA.GT.360.) A = A - 360. 

WD(K-l) = A 


j 
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67 CONTINUE 
GO TO 750 
751 WWE I K- 1 1 =99.9 
W$N?K-lt=99.9 
WSIK-l»=99.9 
WO I K~ 1 ) = 999 .9 
750 CONTINUE 

WIND AT THE SURFACE tS OEFINEO FROM THE BASELINE CARD 

won i = woi 

IFIWD! .60.999. 1 WDI1I*999.9 
WSI U = WSI 

IFIWS! *£0.99.91 WSU)*99.9 

SURFACE HIND COMPONENTS ARE DEFINED 

WWE I l) =WWE I 2) 

WSN( 11 =WS*M 2) 

KT=K T- 1 


PERFORM A 5-Pn(NT WEIGHTED AVERAGF ON 30 SEC «i#IOS-->M ND 
CCMPCNE’ITS ARE SMOOTHED 
TRANSFER STORAGE OF WINOS 


DO 3 25 |=t, KT 
WODI 1 1 = WD( | I 
WSSUMWSUl 
325 CUNT INUft 

RINOmial WEIGHT -COEFFICIENTS ARE USED 

CD 1 1 1 = 0.06 
CD (2 1=0. 25 
COt 3 1=0.38 
CD I 4 1 = 0. 25 
CO t 5 1=0.05 
00 393 1*1, KT 
KT 1 = KT — I 

IF ( KTI . L E. 2 ) GO TO 305 

SMOOTHING DOFS MOT BEGIN UNTIL 2000 M ABOVE GROUND LEVEL 

IF T H INTI II . LT. 2000. 1 GO TO 306 
K 1 = I -2 
K2=l »2 
K4 = 0 
S1=0.0 
S2=0.0 

DO 322 K=K 1 ,K2 
K4=K4*-l 

|F AW "AID VALUE INVOLVED IN THE SMOOTH IS MISSING, DON'T 
SMOOTH THAT POINT 

IF I WWE I K, 1. EC .99.9 1 GO TO 305 
IFIWSNIKI.E0.99.91 GO TO 306 

LOCATE WINDS BASED ON AN ELEVATION ANGLE LESS THEN 9 OG 
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I F t THE r A(K ) .IT. 0.157) IEP(I>=1 
SI = S !*• W’WEI K > *CD I K4 ) 

■ S2=S2+WSN< K )*CD( K4) 

322 CONTINUE 
C 

C INDICATE THAT SMOOTHING WAS PERFORMED (IASM=I> 

C 

I ASM ( I > = l 
WE ( I )=Sl 
KNII )=S2 

IFUBSIS2) .GT.l.OE-5) GO TO 5152 

TT = S 2 

S2=1.0E-5 

IFITT.LT. 0.0) S2=-1.0*S2 
5152 WSU)=SORT ( S 1 + *2*S2**2 ) 

B=Sl/S2 

A= AT AN ( ABS ( B7 >*57. 29578 
IF I B >471 » 4 71 , 472 
4 71 IP(Sl) 473,473,474 

473 A= 36 0. - A 
GO TO 476 

474 A= I 8 0. —A 
GO TO 476 

472 I F ( SI) 473,478,476 
478 A=A+ 180. 

476 CONTINUE 
C 

C MODIFY WIND DIRECTION TO THE GENERAL MET CONVENTION 

C 

A=A*180. 

IFIA.GT.360.) A= A— 36 0. 

WD I T ) s A 
GO TO 393 
306 WD ( I ) = WDDI I ) 

WS( I )=WSS< I > 

WEI I ) = WWE< I > 

WNl I ) =WSN( I > 

393 CONTINUE 


LINEARLY INTERPOLATE WINDS ON BASIS OF TIME TO PLACE ON CONTACTS 
SURFACE VALUES ARE COMPUTED FIRST 


WDDI 1) — WD( 1 ) 

WSSI 1) =WS( 1 > 

IF (WDDI 1 ). £0.999.9) GO TO 838 
IF (WSS I 1 ). EQ.99.9) GO TO 888 
T= 18 0. +KDI 1 ) 

IFtT.GE.36 0.) T= T-3 60. 

Ul 1)=WSC 1) «SINfT/57.29578) 

VI t)=WS(l) *C0S(T/57. 29578) 

GO TO 889 

888 Oil) =99.9 
VI 1 ) =99. 9 

889 R{ 1 ) =0 .0 
ATHI 1) =0.0 
ICT = 1 

DO 613 1=2, 300 

IFITMORI I) .GT.TMNG! KT) ) GO TO 921 
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no 61* J = l C T«300 
K=J 

LOCATE LEVELS ON WHICH TO I NTDRPCLATE 

IFITMNGI J*ll .GF.THOPII II GO TO 615 
61* CONTINUE 
615 ICT = K 

DON'T TRY TO INTERPOLATE MISSING DATA (<MS> 

T FT WOT K) ,E0 .999.91 GO TO 3766 
iriHOI K*l) .EQ.999.9) GO TO 3766 

KEEP UP WITH EL. ANGLES LESS THAN 10 OG AND POINTS THAT HAVE BEEN 
SMOOTHED 

IF I I AS P( K I . EO. I ) I ASMI IM1 
IF T IASMUUI.EO.il I ASMT l I *1 
IF I TfPIK + 1 ) .F0» 1 1 IF PI I 1 = l 
IF(!E D (K).€0.n ifrhmi 

nT = l TMURII ) — T M NG IK I T / 1 THNQ I K* l l-TMNGOO 1 

COMP UT F LOCATION COMPONENTS ON PRESSURE CONTACTS 

XX SI II = XS( K » HU* I XSI K* ll-XSIK) ) 

ZZSI IMZSI KlfilttzslKMI-ISIKI 1 
IF I A ns I ZZS I III .GT. 1.06-51 GO TO 5153 
TT = ZZSI I I 

ZZSI n = 1 .OF-5 * ■ 

in u.lt. o.oi zzsi n=-i.o*zzsm 

COM P'JtE RANGE 1 Rl ANQ AZIMUTH I ATHI OF SONDE 

5151 IFtXXSI U.NE.ZZSm I GO TO 3033 
ITIXXSM ). EO. 99.91 R 1 1 I* 999. 9 
IF (XXSI I |. 60.99.91 ATM III *999. 

GO TO 303* 

303 3 Rill si SORT | I XXSI I) )**2*U2S< IT I **211/1000. 

F=XXSI ll/ZZSI I I 
0= AT AN ( ABS | F 1 1 *57. 29 570 
IFTFT67L1, 6711,6721 
67L1 IFtXXSmi 6731,6731.67*1 
6731 0=360.-0 

GQ TO 6761 
67*1 0=180.-0 

GO TO 6761 

6721 IFIXXSMJ1 6781,6781,6761 
6781 D=0*L80. 

6761 0=0*180. 

1FI0.GT.360.) 0=0-360. 

ATH( I MO 
303* CONTINUE 
C 

C COMPUTE NEW WIND COMPONENTS, DIRECTION ANO SPEED 

C 

U( I MWE< K) +DT*< WEI K*ll-WE(KI 1 
VlII-WNIK) ♦DT*1WN(K*1I-WNIKM 
IF (ABS I VC I II. GT. 1.06-51 GO TO 51 5* 

TT = V 1 1 I 
VI IM1.0E-5 
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371 
3 73 

3 74 

372 
376 
3 76 


ifitt.lt. o.oi v(n=-i. 0 '»vm 

I LSE WSS AND HOD FOR SPEED AND DIRECTION NOW 

5154 V»SS( 11 =SQRT( (U( H I **2T.(Vm )**2) 

R=U( I) /V< I ) 

A=ATAN( ABSt B 1 1 * 57 . 29 57 8 
IF1B1 371,371,372 
IFIUUM 373,373,374 
A- 36 0. - A 
CO TO 376 
A= 1 8 0. - A 
GO TO 376 

IFIU1I1I 378,378,376 
A=AH80. 

A= A *-18 0 . 

IP! A.GT.360. 1 A = A- 360. 

HDD! I) =A 

MODIFY WIND OIRECTION TO THE GENERAL MET CONVENTION 
wntll I) =WDD< I 1 + 180. 

IF1WD01 I 1.GE.360.) MOD 1 1) =WD0( 1 1 -3 60 . 

MAKE U POSITIVE POINTING EAST AND V POSITIVE POINTING NORTHWARD 

■ U( 1 1 =- 1. 0*U I 1 1 • 

V(I)=-L.O*VU> 

GO TO 613 
3766 W00( I) =999 .9 
WSS( 11=99.9 
U( 11=99.9 
VI 11=99,9 
IF P < 11 = 0 
Rt I 1=999.9 
ATH1 11=999. 

613 CONTINUE 


PRINT THE RESULTS 

921 KH = 0 

IFIKAT.LT. 10) KT =K AT 
00 5 5 K=l, KLM 

ADO STATION ELEVATION TO HEIGHT CF SONDE 
YS IK )=YS(K 1 +HT 

PRINT 9’S IF WIND VALUES END BEFORE THERMO DATA 


IF(TMORIK) 
IF ITMOF ( K) 
IF < THOR IK) 
IF (THOR I K) 
IF ITMQP 1 K) 
IF (TMORI K) 
IFITMOR(K) 
IF IT HOP. ( K) 


GE.TMNGIKT- 
GE.TMNG ( KT- 
GE.TMNGl KT- 
GB.TMNGIKT' 
GE .TMNG ( KT- 
GE.TMNG1KT- 
GE .TMNG 1 KT- 
GE.TMNG I KT- 


1)1 U 1 K) =99. 9 
1)1 VI K) =9 9.9 
'll 1 R ( K) =999.9 
1)1 AT HI K) =999 . 
Ill WDD1 K 1 =999 .9 
1)1 WSS1 K) =99.9 
111 IFP( K> =0 
111 I A SM ( K ) = 0 


1 
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C . PRINT A * T C SHOW SMOOTHING AND INTERPOLATION INSTEAD OF A 1 
C 

IT II OR IN IK I . EO . 0 ) lORf NIK) =B LANK 
|Ff IOR IN1KI .EO.tl IORIN(K|=AST 
irHFP.lkl.EQ.OI IFPIK)=BLANK 
IF | | FPIK). EQ.U I F P I K ) = A S T 
IM I ASM! K> . EQ.OI I A SMI Kl = BLANK 
I r | ! A c M IK), ED • 1 ) I A SM| K I =A ST 
IF (MODI KHt A 5) I 9A,91,9A 

93 WR ! T El 6,90 1 ID 5, I NAME! 1 1 , 1 =1 ,7 ) , II) l , IX, I Y . 1 1 , 1 03 . 1 OA ,KLM,P I KLM ), 

1 IlMtN 

90 FORM AT ( 1H1 , 56X, 'STATION NO. • , I 5 , / , 5 1 X , T AA , / / , 

I 5 7X, I2,2X, 3 A3, ' 19' , I 2, /, 61 X, I A ,1 X, 'GMT* ,50X, 13 , 1 X , F5 .0 , 2X, III 
IF ( JlMlN.NE .0) MRI TF (6,9771 

977 FOR” ATI ' *, 'ANGLES ON THE HALF MINUTE HAVE BEEN LINEARlV INTERPOLA 
ITFt) FROM WHOLE MINUTE VALUES'! 
f'P IT El 6,96 21 

9 62 FORMAT I 'O' , • TIME ' , AX , ' CN TCT ' , A X, 'HEIGHT' , AX, 'PRES', 

1 5X ,' TEMP' ,3X ,' DEW P T ' , AX » ' 0 1 R ' » AX, 'SPEED', 3X, 

2'U COMP ' ,3 X , 'V COMP* ,3X, 'POT T',3X,»E POT T'» 

3 3;,, 'MX H TO • , SX.'RH', AX, 'RANGE' ,2X, 'AI'» 

V,o IT El 6, 96 3 ) 

963 FORMAT!' • , IX, 'MIN* , 15X, 'GPM ' , 6X,'MB',6X,'DG C • , 

LAX, ' DG r • , 6X, * DG' , AX, ' M/SEC* ,AX, 'M/SEC*, AX, 

2 'M/SEC ' ,AX, *0G K'.AX.'DG K ' , 6X , ' GM/KG' ,AX, 

3 1 P C T ' , 6X , ' KM ' , 3X , ' DG ' , / ) 

9 A WPITEI6.60I TMORIKI ,CTC(K» , V SI Kl ,P(K),TCIK), IORtNIKl , 

1 TP IK ) , W3DI K ) , WSS IK ) , IF PI K I ,U IK I, VI K) , PTK I K > , EPOT IK » , 

2 W ( K I ,H M ( R( , P I K ) , ATH I K) 

60 FPP M AT I ' ' , F9. 1 , 3X,F5. I, 3X,F7. |, 3X,F6. 1 ,3X*f'5. I» 

1 AA , F 5. 1 , 3X . F 5. 1 , 3X, F 5. 1 , A 1 , 2 X, E6 . 1 , 3 X, F6 . 1 , 

2 3X.F5. 1, AX, F5. 1, 5X.FA. I, AX.F5. 1.3X.F5. 1 , 1X.FA.0I 
KH*KHf"l 

55 CONTINUE 
C 

C FP ITE THE ! AFFIRMATION ON TAPE USING NO FORMAT 

C 

FR IT El 3) 105, NAME, 101, IX, IV, 12,103, IDA, 

1 KL M, P( KL M| , I l MI N, TMOR, CTC, YS ,P ,TC, 

2 IP® I N, ID, HDD ,WSS,IFP ,U , V , PTK, EPOT, 

3 W, HM ,R , ATH 
999 CONTINUE 

GO TO 9266 
9366 WP 1 T El 6, 9367) 

9 367 FORMAT! • 1' , 'ERROR IN READING TAPE'! 

9266 CONTINUE 
REWIND! 

REWIN03 

STOP 

END 

C 

c 

SUPP. CUT I NE MOPRTIIX, IA, I B, IC I 
DIMENSION I Mil 12), (M2I 121, 1M3I 12) 

CAT A I M l /3HJAN,3HFEB» 3 MM AR , 3HA PR , 3HMA Y , 3HJUN , 3HJUL , 3HAUG.3HSEP , 

I 3HQCTY3MN0V ,3HDECA 

DATA I M2/3 HUAR , 3HRUA ,3HCH ,3HIL ,3H ,3HE ,3HY ,3HUST, 
l 3HTEM, 3H0B E, 3HEMB, 3HEMB/ 

CATA I M3/3HY , 3HR Y ,6»3H , 3HBER ,3 HR ,3HER , 3HER / 

T I - I X 

IA = tMtlin 
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IB= I M2 ( I I) 

I,C = TM3 tIU 

RETURN 

END 

C 

C II ID AT A - 

C 

54 

11001 Cl 52 3466 08670 

22001 0362 0524 09746 

22002 C362 3466 09831 

22003 C449 3479 09769 

22004 0423 3515 09847 

22005 C4 5 1 3510 09796 
20 L 0003 2458 08168 
202 C004 2583 08027 
208 CO i 3 3291 08002 
211 0008 2792 08252 
213' 0044 3125 0S241 
221 CO 22 3050 08652 
226 CO 5 7 3225 08641 
232 C001 2933*06941 
235 OICO 3233 09008 
240 CO 05 3 00 8 09318 
248 CG 79 3250 09382 
250 C0C7 2592 09743 
255 C033 2833 09692 
259 0399 3233 09825 
261 03 14 2933 1C092 
265 C 8 7 3 3 193 10219 
304 C004 0526 07537 
311 0246 3394 08332 
317 C27 5 3608 07995 
327 Cl 80 3625 0,3657 
340 0079 3474 0V224 
349 C438 3690 09591 
363 10 95 3524 IC1C-9 
402 0004 3784 07549 
405 0065 3898 07748 
425 C246 3836 08254 
429 0298 3986 08410 
433 G 1 7 5 3866 08898 
451 C 7 91 3 776 09998 
456 C263 3907 09562 
486 0007 4066 07377 
494 0016 4167 06998 
518 C086 4275 07378 
520 C359 4053 08024 
528 C21S 4293 07874 
532 C2 00 4067 08968 
553 C4 C 3 4136 09601 
562 C847 4114 1 0068 
606 0020 4366 07032 
637 C23o 4297 08375 
645 0210 4449 08815 

654 C3 92 4439 09820 

655 0316 4558 09418 
66? C966 4405 1 0307 
712 C 1 9 1 4686 06801 
734 C22 1 4686 07435 
76:4 C5C3 4676 10075 

747 C359 4857 09337 

C END 


MARSHALL SPACE FLIGHT CENTER 
NORMAN, OKLA 
FT. SILL, OKLA 
LINDSAY, OKLA 
FT. CD 88 , OKLA 
CHICKASHA, OKLA 
KEY WEST, FLA 
MIAMI, FLA 
CHAP. LF ST ON, SC 
TAMPA, FLA 
WAYCROSS, GA 
EGLIN AFB, FLA 
MONTGOMERY, ALA 
BOOT HV I L LE , LA 
JACKSON, MISS 
LAKE CHARi.ES, LA 
SHREVEPORT, LA 
BROWNSVILLE, TEX 
VICTORIA, TEX 
STEPHENVILLE, TEX 
DEL RIO, TEX 
MIDLAND, TEX 
HATTER AS, NC 
ATHENS, GA 
GREENSBORO, NC 
NASHVILLE, TENN 
LITTLE ROCK, ARK 
MONET TE, MQ 
AMARILLO, TEX 
WALLOPS ISLAND, VA 
DULLES AIRPORT, VA 
HUNTINGTON, WVA 
DAYTON, OHIO 
SALEM, ILL 
DODGE CITY, KAN 
TOPEKA, KAN 
KENNEDY AIRPORT, N Y 
CHAT AM, MASS 
AL8ANY, N Y 
PITTSBURG, PA 
BUFFALO, N Y 
PEORi A, ILL 
OMAHA, NEB 
NORTH PLATTE, NEB 
PORTLAND, ME 
FLINT, MICH 
GREEN BAY, Wl$ 

HURON, S D 
ST CLOUD, MINN 
RAPID CITY,* S D 
CARIBOU, ME 
SAULT STE MARIE, MICH 
BISMARCK, N D 

INTERNATIONAL FALLS, MINN 


4 

1 
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APPENDIX E 

Use of the First Difference Program 

The number of soundings beginning with the first on which first 
differences are computed is read from a card; results from the master 
reduction program are read from magnetic tape using DCB = (REGFM = VS, 
DEN = 3). The program has an object code of 8,240 bytes and an array 
area of 33, 148 bytes. Results are printed using the same format that 
was used for the master reduction program (Fig. 6). 
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THIS PROGRAM COMPUTES FIRST DIFFERENCES OF 
QUANTITIES COMPUTED AT PRESSURE CONTACTS 
SUMMER 1 97 4 — HENRY FUELBERG 
ARRAYS A AND IA ARE ORIGINAL SOUNOINf DATA 
ARRAYS B AND IB A? E FIRST DIFFERENCE DATA 

DIMENSION A (230, 16) ,B( 230,16), I A (230, 21 , IB(230,2) 

INTEGER AST, BLANK, NAME (7) 

DATA AST, BLANK/ 1H*,1H / 

RE HINDI 

N IS THE NUMBER OF SOUNDINGS TO COMPUTE 
RFAO, N 

DO 500 I SOUND= 1 , N 
READ IN THE DATA 

RF AD ( 1 , ERR = 9 l, END= 50 II f D5, N AME , ID 1, IX, I V, I Z , I D3 , IDA , 
1,KLM, PKIN.I l MIN, ( <A( I ,j ) , ! = i, 23 0) ,J=1,5) , <IA< I, 1) ,1*1,230), 
2 ( ( A( I, J) ,1 = 1,230), J = 6, 8) , f I A II .21,1 = 1,230), 

3 t t A ( I, J), 1 = 1,230), J = 9, 16) 

COMPUTE FIRST DIFFERENCES 

klmi=klm-i 

COMPUTE 8Y ROWS 

CO AGO 1=2 , KLMl 

INITIALIZE IE ARRAYS WITH BLANKS 

I B ( T ,1 ) = BL A NK 
I B ( t ,2 ) = BL A NK 
K= I 

COMP UT F BY COLUMNS 

00 300 L =1 , 16 -V 

DON'T COMOUTE FIRST DIFFERENCES IF A VALUE WAS MISSING (9'S) 

IF(A(K,L).EC.999.9)G0 TQ 313' 

IF< A(K,Li.F C.99.91G0 TO 314 
IF (AIK-1, L) .E0.999.91G0 TO 313 
IF(A(K-1,L) .E0.99.91GO TQ 314 
B ( K,L)=A(K,L>-A(K-l»L) 

GO TQ 300 1 1 

313 B(K,L) = 999.9 
GO TO 300 

314 8(K,L)=99. 9 i 

30 0 CONTINUE - ■ 


KEEP UP WITH STARRED VALUES 
INDICATE INTERPOLATED VALUES BY A * 

INDICATE WINDS; THAT WERE COMPUTED WITH ELEVATION ANGLES LESS THAN 10 DG 

IF ( I A( K» 1) . EQ. A ST) I B IK , 1 ) = AST 
IF ( IA( K-l, 1 ) .E Q. AST) IB (K, l)=AST 
IF (I A(K, 2) .EQ.AST) I B(K,2)=AST 
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IF! I A(K-i, 2I.EQ.ASTI IB (K, 2 CAST 

<4oo continue 

C PRINT THE RESULTS AS BEFORE 
KH=0 

00 55 K=2,KLM1 
IF IMODIKH, ASH 94, 93,94 

9 3 HR IT El 6. 90) IDS, I NAME II I, 1 = 1.7) , 101 , IX, I Y, I Z, 1 03,. 1 04, 

IKIH, PMIN.I 1M1N 

90 FORM AT 1 1 HI , ' **F I RS T 01 FFERENCE St*' ,36X, • STAT ION NO. ', 
115,/ ,51X,7A4,//,57X. I2,ZX-,3A3, *19' ,1 2,/,6tX, 14, IX, 

2 'GMT SOX, I3,1X,F5. C,2X, 1 1 > 

IF I I IM I'il •NE.0)WRITF(6,977) 

977 FORMAT ( •, 'HALF MINUTE ANGIES WERE LINEARLY INTERPOLATED*! 

WR I T E ( 6, 96 2 I 

962 FORMAT t *0' , * ?! ME • , AX ,* CN TCT* ,4X, •HEIGHT* ,4X, 'PRES*, 

15X. ' TEMP* , 3X.’ DEW P 1 • * AX , • OT R* ,6 X, • S PEED • , 3X , 

2'U C CM P ' , 3 X , *9 COMP* ,3X, 'POT T',3X,'E POT T' , 

33Xf'MX RTO • ,5X, • RH* , AX ,' RANGE' »2X, *AZ' I 
WRITE 16,9631 

963 FORMAT!' • , IX » 'M IN ' * 15 X * 'OPM * * 6X * ' MB ' «6X » • DG C • , 

1 AX , ' DG C ' t 6 X » • OG • * AX » ' M/SEC- ' ,4 X, *M/SEC * , 4X* 

2 'M./S EC ' ,4X , *0G K*,4X,*DG K ' , 6X , • GM/KG' , AX, 

3 'PCT ' * 6X,* KM* , 3X , 1 DG ' , / ) 

C 

C THE PRESSURE CONTACT AIK, 2) IS NOT DIFFERENCED BUT THE FORWARD 
C VALUE IS PRINTEC IN" ORIGINAL FORM 
C 

9A WRITFI6.60) 8(K,1),A(K,2),(B(K,L),L=3,5I,IB(K,1), 
t IRlK,L),t=6,ai,IBIK,25 , I B i K, LI , L=9 , 16) 

60 FO p MATi * * ,R5. 1,3X,F5. l,3X,F7. 1, 3X,F 6.1 ,.3.X, F5. 1 , 
t AA.F5. 1,ZX,F6.1,3X,F5. I , A1 ,2 X, F6 . 1 , 3 X, F6 . 1 , 

23X.F 5. 1, AX, F5. I, 5X ,F A. 1 , AX , F 5. I , 3X ,F 5. I , IX,F A. 01 
KH*K Hf 1 
55 CONT INUE 
500 CONT INUE 
GO TO 501 
91 HR 1 T El 6, 91 1 I 

911 FORMAT, • 1* , 'ERROR IN READING THE TAPE*! 

GO TO 501 
501 REWIN01 
STOP 
END 
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APPENDIX F 

Use of the Correction Program 

Two tapes are needed; the first contains the original angle, 
ordinate, and baseline data. The tape is non-labeled and has DCB = 
(RECFM = FB, LRECL = 80, BLKSIZE = 1600, DEN = 3). The second tape 
will contain the corrected data and has the same tape parameters as 
the first. The program has an object code of 7,040 bytes and an array 
area of 40,000 bytes. The first data card read specifies the total 
number of soundings for which data are to be transferred to the new 
tape. This card is followed by a sequence of cards for each sounding. 
One card specifies the number of corrections to be made on each sounding 
which is followed by a series of locator cards and correction cards that 
were described in Section II-E. If no corrections are made to a 
sounding, only the date and time of the sounding is printed on paper, 
but if corrections are made, the new sounding data is printed until two 
records after the final correction. The printed output facilitates a 
verification that the corrections have been made. 
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TH t s PROGRAM replaces incorrect records on the rah data tare 
WITH THE correct records 

SUMMER 1 97 A — HENRY FUBLBERG 

CHAP AC TER* 80 OCARD I 500 I , NCARD 
REWIND 1 
REWINDS 

NSnUND IS THE NUMBER OF SOUNDINGS TO CHECK 

REAO »NSOUND 
00 100 I SDUNO= 1 i NSOUND 
C 

C NCOR IS THE NUMBER OF CORRECTIONS IN THE GIVEN SOUNDING , 

t 

READ, NCOR 
C 

C P E AO THF * LEADER CARD* 

C 

READ 1 1 . 102 , EPR = 90i> » IMO, IDA , I VR, I TM, I STA , ICDS 
102 FORMAT ( 12, I X, I 2 , IX, I 2, IX , I 4, IX, 15, 2X, 13) 

C • 

C IF CORRECTIONS ARE TO BE MADE ON THE 'LEADER CARD', 

C FERE IS THE PLACE TO DO IT 

C 

WRITE! 3, 10 2 1 1M0,IDA,IVR,ITM,I STA, ICDS 

C 

C WRITE THE 'LEADER CARO' ON TAPE ANO PAPER 

C IF NO CORRECTIONS ARE TO BE HADE ON THIS SOUNDING, 

C THIS IS THE ONLY PRINTED INFORMATION THAT WILL APPEAR 

C FOP- THIS SOUNDING 
C 

WRIT El 6, 10 21 IMO, IDA, I YR , I TM , I STA, ICDS 
951 CONT INUF 

C - - 

C ICDS IS THE NUMBER OF RECORDS IN THE SOUNOING 

C READ ALL THE RECORDS IN THE SOUNOING 

C 

00 200 1=1, ICDS 
C 

C OCARD IS THE RECORD CURRENTLY ON TAPE 

C ■ ■ ■ / 

RE AD ( l , 202 < ERR = 906 1 OCARD! 1 1 


oooono oooo nnonooooo r> r> r> 
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20 2 FORM AT (480) 

20 0 CONTINUE 

IF < NCOR. NE . 0 ) GO TO 300 

WRITE THE RECnPDS ON A NEW TAPE IF NO CORRECTIONS ARE NEEDED 

CO 250 1 = 1, ICDS 

HP IT El 3, 25 1 ) OCAROl 1 ) 

251 FORMAT (A 80) 

250 CONTINUE 
GO TO 100 
300 1ST" 1 

MAKE SUBSTITUTIONS FOR INCORRECT RECORDS 
CO 400 J= l » NCOR 

CAPONO IS THE NUM3ER OF THE RECORD TO CORRECT 

reaDjC ardno 

NCARD IS THE CORRECT RECORD TQ BE PLACED ON TAPE 

PE AD 1 5 » 4 01 ) NCARD 
401 FORMATlASO) 

CO 450 K=t ST, ICDS 
KL=K 

IF (K.EQ.CARDNO) OCARDI K) =NCA RD 

WRITE ON TAPE AND PAPER ALL RECORDS UP TO THE LAST CORRECTION 
OF THE SOUNDING 

KPITE(6,48l) OCARD(K) 

481 FORMAT ( l X.A 80) 

WR I T E ( 3* 45 1 ) OCARD(K) 

451 FORM AT (A 80) 

IF (K.EQ.CARDNO) GO TO 455 
450 CONTINUE 
455 CONTINUE 
I ST = KL ♦ l 
400 CONTINUE 

IF C I ST.GT. ICDS) GO TO 100 

WRITE THE REMAINING RECORDS ON TAPE 

WR ITE( 3» 45 1 > (OCARU(KZ) ,KZ=I ST, ICDS) 

WRITE TWO MORE RECORDS ON PAPER AFTER THE LAST INCORRECT RECORD 

WP ITE ( 6,48 1 ) OCARDI I ST) ■ 

WRITE! 6 » 48 1) OCAROl IST + L) 

100 CONTINUE 
GO TO 001 

906 WR I TE( 6,907 ) ' " ■■ ■, ' 

907 FORMAT ( ’ l* , ’ERROR IN READING TAPE’) 

901 REWIND1 

REWIND3 : 

STOP 

FND 
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APPENDIX G 

Use of the Pressure Interpolation Program 

Data from the master reduction program is read from a non-labeled 
tape with DCB = (RECFM = VS, DEN =3). A card is read that specifies 
the number of soundings that will be interpolated beginning with the 
first sounding on tape. The program has an object code of 9,768 bytes 
and an array area of 20,188 bytes. Results are printed on paper with 
the same format that was used in the master reduction program (Fig. 8), 
and they are transferred to a second non-labeled tape with DCB = (RECFM 
VS, DEN =3). The format of the tape is described in Table 7; quantitie 
have a dimension of 50 instead of 230 which was used in the master 
reduction program. 
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CONVERT DATA TO 25-MB INCREMENTS 
SUMMER 197 A— HENRY FUELBERG 


B IS THE 2 5-MB DATA 
A IS OR IGI ONAL DATA 

'.I ME NS ION A ( 230, 16) . B(50,16>, IAI230,2>, IB(50,2> 
INTEGER AST, BLANK, NAME! 71 
CATA AST, BLANK/ IH*, IH / 

REWIND! 

REWIND3 

N IS THE NUMBER OF SOUNDINGS TO COMPUTE 
READ.N 

00 500 LS= l,N 

INITIALIZE THE NEW B ARRAY 

DO 5 03 1=1,50 

00 509 J=1 , 16 
B( I, JJ*0.0 

509 CONTINUE 

1 B ( I ,1 ) = BLANK 
I B ( I ,2 > = 8L ANK 

503 CONTINUE 


READ IN THE DATA 

READ 1 1 , FRR = 9 1 , ENO= 501 ) I D5.NAME, ID 1 , IX , I Y, I Z , I D3 , 1 09 , 

1KLM, PM IN, I 1MIN, l (A( I , J ), 1 = 1,230) ,J = l ,5), ( I A< It II ,1 = 1,230), 

2 II A< I, J > ,1 = 1,230) , J = fe, 8) ,( IAU ,2), 1 = 1,230) , 

3 (( A( I, J > , I = 1,230), J=9, 16) 

C 

C LEVELl ON NEW SCHEME I S THE SU.RFACE 

C 

CO 200 1 A= l, 16 
B(1,LA)=AC l.LA) 

200 CONTINUE 

181 1 ,1 ) = IA( 1,1) 

IB ( l , 2 ) = I A (1,2) 

C 

C BEGIN LINEAR INTERPOLATION STARTING AT 1000 MB 

C 

P=1000. 
no 300 1=2,91 

c 

C IF SURFACE PRESSURE IS LFSS THAN THE PARTICULAR 25-MB LEVEL, 

C PRINT 9 • S INSTEAD OF VALUES 

C 

IFIP.GF.Al 1,9)) GO TO 550 
C 

C IF THE 25-MB PRESSURE IS LOWER THAN THE MINIMUM SOUNDING PRESSURE, 

C PRINT 9 1 S INSTEAD OF ACTUAL VALUES 

C 

IF(P.LT.PMIN) GO TO 550 
ICT= 1 
C 

C FIND PRESSURE LEVELS ON WHICH TO INTERPOLATE 
C 

DO 301 J=I CT.230 
K=J 

IF(A(J + 1,9 ) -LE.P) GO TO 302 
301 CONTINUE 


302 1C T * K , 

B< !, AMP 
C 

C DP I S THE INTERPOLATION FRACTION 

C 

DP = ( P-A< K, A ) )/( A IK* l, A l-AIK, All 
CO 305 I -t . 16 
C 

C DON'T INTERPOLATE PRESSURE IL.EQ.AI 

C 

IF (L .EQ. Al GO TO 305 
C 

C IF A VALUE USED IN THE INTERPOLATION IS HISSING PRINT 9'S 

C 

IFIA(K*t,L).FQ.999.9) GO TO 313 
IF (AtKU'Ll .EQ.99.91 GO TO 3 1A 
| F ( A I K. , I _I.FC.999.9I GO TO 313 
It (AIK. LI. EC. 99.91 GO TO 31A 
C 

C A SPECIAL SCHEME IS USED TO INTERPOLATE WIND DIRECTION 
C IL.EO.T) SINCE IT MAY OSCILLATE AROUND 360 DEG 

C • 

IFIL.EO.T) GO TO 320 

IF< ( L.F0.L6 I .AND. I A(K* 1,161. EO. 959,11 GO TO 326 
IF((L.FC.16).AND.(AIK, L6I.E0.999.il GO TO 326 
Bl I. U» “A IK « L I ♦IJP*I A IKE 1, l.l-A IK,L II 
GO TO 305 
313 Bt I , L) =999 . 9 
GO TO 305 
3 1 A 8( I,LM99.9 
GO TO 305 
32A Eli, 161 = 999. 

GO TO 305 

320 CAZ=AIKH,LI-AIK,L) 

E ( I » LMA (K , L )>OP*OAZ 

IF I A BS I DAI » . GUI 00. I B 1 1 , L » = B( I , L) H BO. 

IE I B ( I .U. GE.3 60.1 BIT » L MB! I, LI-360. 

305 CONT INUF 
C 

C KEEP UP WITH STARRED OUANTITIES 

C 

IF ( I At KM, LI .EO.AST) IBII , II ’AST 
IF ( 1 AIK, II. EO. ASTI t B( I.IMAST 
IF! I A( Kf I, 21 .FO.AST) I B( I ,2) *AST 
IF ( I AIK, 21 .EO. ASTI I B( 1 , 2 I =A ST 
GO TO 309 
C 

C FILL IN MISSING DATA WITH 9*5 

C 

550 CO 307 L-l.13 
B I I , L ) = 99. 9 
30T CONT INUE 
B I I, AMP 
B ( I , 12 M 999.9 
Bt I, LA 1 = 999.9 
Bl 1 , 15 1 = 999.9 
811,161=999. 

C 

C SUBTRACT 25 MB FROM P AND INTERPOLATE AGAIN 
C 

309 P*P-25. 
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300 CONTINUE 
'C 

C PRINT THE RESULTS 

C 

WR IT El 6, 90 > 105, (NAME! I) ,1 = 1,71, 101, IX, I Y, I Z , I 03 , 1 DA ,KLM ,PM IN» 

1 1 1HIN 

90 FOPMATUHl , 56X, • STATION NO. • , I 5 , / ,5 l X , 7 AA , / / , 

l 57X, !2,2X, 3A3, '19' , !2, / ,6LX, IA, LX, 'GMT' ,50X, 1 3 , 1 X, F5 .0 , 2X, 1 1) 

IF! I 1MIN.NF .0) WRITE (6,977) 

977 FORMAT! • 'ANGLES ON THE HALF MINUTE HAVE BEEN LINEARLY INTERPOLA 
l TEO FROM WHOLE MINUTE VALUES') 

WRITE! 6, 96 2) 

962 FORM AT ( 'O' , • T I ME • , AX , ' CN TC T * , A X , • HE I GHT ' , AX , • PRE S • , 

1 5X , ' TEMP • , 3X , * DEW P T ' , AX . ' 01 R • , AX, ' SPEED' ,3X, 

2 'll COM P ' , 3 X * ' V COMP' ,3X, 'POT T',3X,'E POT T* , 

3 3X. ' MX R TO • , 5X , ' RH 1 , AX , ' RANGE* ,2X, ' AZ’ ) 

WR ITE! 6,963 ) 

963 FORMAT!' ' , IX, ' K IN ' , 15 X, ' GPM • , 6X , • MB • , 6X , ' DG C ' , 

1 AX,' OG C',6X, 'OG'.AX,' M/SEC' ,AX, 'M/SEC', AX, 

2 ' M/S EC ' , AX , 'OG K'.AX.'DG K« , 6X , ‘ GM/KG* , AX. 

3 ' P C T ' , 6 X , ' K M * , 3 X , • D G ' » / J 
DO 55 IZ=L,A1 

V.P ITE! 6,60 > ( Bl IZ» I ) , I =1, 5), IB (I Z, 1) ,<B( IZ, I J, 1 = 6,8) , 
l IB (l Z, 2) , ( B t IZ, 1 ), 1 = 9, 16) 

60 FORMAT!' ' , F5. 1.3X.F5. 1,3X,F7. 1 , 3X , F 6 . 1 , 3X , F 5. 1, 

1 AA,F5.1,3X,F5.1,3X,F5. 1 , A l , 2 X , F6 .1 , 3X, F6 . 1 , 

2 3X.F5. l,AX,F5.l,5X,FA. 1,AX,F5. l,3X,F5.ir IX.FA.O) 

55 CONT INUE 

WR ITE! 3) t05,NAM&, ID 1, I X , I Y, I Z , I D3 , 1 DA , 

1KLM, PMTN, I IMIN, ( ID! I ,J ) , I = L, 050) ,J = 1 ,5) , (I 3! I, 1) ,1 =1,050), 

2 ! t B( I, J) ,1 = 1,0 50), J = 6, 8) , (IB (I ,2), 1=1,050) , 

3 ( ( 3! I, J) ,1 = 1,050), J = 9, 16) 

500 CONTINUE 
GO TO 501 

91 WRITE(6,9ll> 

9 11 FORMAT .('l' , 'ERROR IN REAOING THE TAPE') 

501 RFWIND1 
REWIND3 
STOP 
END 


i ' ' ' » 


I 
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APPENDIX H 

Known Errors Remaining in the Reduced Data 


Locations of errors in the basic parameters — temperature, 
pressure, humidity, time — are given by date and time. Errors in the 
basic parameters will be reflected in the qtiantities derived from them 
such as pressure-altitude, winds, potential temperature, etc. Suggested 
corrections are given where possible for some of the derived parameters; 
other corrections can be calculated from the corrected basic parameters 
and the appropriate equations. 


Station Date /GMT 

221 All time 

Eglin AFB, periods 

Florida 


250 12/0600 

Brownsville, 

Texas 


Azimuth angles are 180° out of phase. 
Correct derived wind direction and 
balloon azimuth location by 180°. U and 
V wind components are 180° out of phase. 

The baseline (surface) wind direction 
should be 140°. Correct U and V wind 
components accordingly . 


260 

Stephenville, 

Texas 

All time 
periods 

SEP on the raw data tape is indicated as 
station 259 instead of station 260, The 
error does not exist in other tapes. 

261 

Del Rio, 
Texas 

11/1500 

The surface pressure should be 966.9 mb. 
Pressure-altitude may be corrected by 
subtracting 268 m from each value given. 

261 

Del Rio, 

11/2100 

The surface wind direction should be 330°. 
Correct U and V wind components accordingly 


Texas 

The surface pressure should be 1013.7 mb. 
Pressure-altitude may be corrected by- 
subtracting 34 m from each value given. 


494 12/1200 

Chatam, 

Massachusetts 
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Station 

Date/GMT 


520 

Pittsburg, 

Pennsylvania 

11/1800 

The surface pressure should be 968.8 mb 
at contact 8.2. Correct pressure-altitude 
by subtracting 104 m from each value given 
Contact 8 is non-existent. 

520 

Pittsburg, 

Pennsylvania 

12/1200 

The surface pressure should be 961.3 mb. 
Pressure-altitude may be corrected by 
subtracting 21 m from each value given. 

528 

Buffalo, 
New York 

12/0900 

Abrupt change in elevation angle at 46 
min after release. Cause unknown. 

637 

Flint, 

Michigan 

11/1500 

The surface pressure should be 979.3 mb. 
Add 52 m to correct pressure-altitude-. 

734 

Sault St. 

Marie , 
Michigan 
747 

International 
Falls, Montana 

All time 
periods 

All time 
periods 

Sondes were released during light rain 
and/or fog in near freezing temperatures. 
Very high humidity values may be due to 
a faulty sensor, and cannot be corrected. 

11001 
Marshall 
Space Flight 
Center, Ala. 

All time 
periods 

Incorrect station elevation was used; 
subtract 12 m from all heights. 

22004 
Ft. Cobb, 
Oklahoma 

12/0100 

The surface pressure should be 961.7 mb. 
Add 93 m to all heights to correct 
pressure altitude. 

All 

stations 

All time 
periods 

Station number, date, and time should not 
be read from ordinate or angle records on 


the raw data tape since these values are 
often incorrect. These values should only 
be read from the baseline record. 


| 

[ 
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APPENDIX I 
Missing Soundings 

Soundings were not computed at the following stations and times for the 
stated reasons. Sounding;? are available at other stations for each of 
the 9 -time periods. 


Station 

Date /Time 

Reason for Omission 

208, Charleston 

12/0252 

Technical problems in the reduction 

process. 

226, Montgomery 

11/1500 

Ordinate data not available due to 
a malfunction in equipment. 

255, Victoria 

12/111,5 

Technical problems in the reduction 
process . 

265, Midland 

11/1200 

Ordinate data not available due to 
a malfunction in equipment. 

22003, Lindsay 

12/0300- 

12/1200 

Soundings not taken. 

22004, Ft. Cobb 

12/0300- 

12/1200 

Soundings not taken. 

22005, Chickasha 

12/0300- 

12/1200 

Soundings not taken. 
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